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ABSTRACT 


The  effect  of  earth  cover  on  the  far-field  fragment  density  expected 
from  the  accidental  detonation  of  stored  munitions  was  estimated  by 
preparing  three  models  of  fragment- cover  interaction.  Comparisons  of  the 
theoretical  calculations  with  limited  experimental  data  show  that  the  model 
wherein  the  crown  of  the  earth  cover  does  not  retard  any  fragments  gives 
the  best  agreement.  Models  for  fragnent-fragment  interaction  which  effec¬ 
tively  account  for  stack  configuration  lead  to  a  simplified  model  for  the 
effective  number  of  munitions  contributing  to  the  far-field  fragment 
density.  An  approximation  technique  for  the  rapid  calculation  of  the  far- 
field  fragment  density  was  prepared  to  assist  in  the  ready  evaluation  of 
any  model.  Tentative  quantity-distance  relationships  for  four  munitions 
were  prepared.  Parametric  studies  of  the  effect  of  altered  mass  distribu¬ 
tions  and  fragment  shape  were  conducted  to  assess  possible  differences 
between  accidental  detonation  source  parameters  and  arena  data  source 
parameters. 
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Chapter  I 
INTRODUCTION 


1.1  Objectives  of  This  Study 

A  theoretical  and  analytical  examination  of  the  far-field  fragment 
distributions  from  the  accidental  detonation  of  stored  munitions  was  under¬ 
taken  to  define  and  illustrate  the  dependence  of  the  fragment  hazard  on  the 
many  parameters  that  describe  the  munition  store.  High  explosive  bombs 
and  projectiles  may  be  stored  in  rectangular  block  stacks  containing  a 
wide  range  of  explosive  quantity,  between  earth  mound  barricades  in  the 
open,  in  above  ground  magazine  structures,  or  in  earth-covered  igloos. 
Explosive  quantities  up  to  250,000  lbs  are  permitted  in  above  ground 
stores,  and  up  to  500,000  lb  in  earth-covered  magazines.  The  present 
design  of  the  revetment  or  igloo,  as  far  as  fragment  retardation,  prevents 
the  fragments  from  striking  an  adjacent  munition  store  and  detonating  it 
by  impact.  This  study  is  concerned  with  the  fragment  hazards  beyond,  say, 
one  thousand  feet  from  the  munition  store  and  concentrates  on  those  frag¬ 
ments  which  are  dangerous  to  personnel,  vehicles  and  structures  there. 

The  retardation  effect  of  the  igloo  is  the  primary  interest  of  this 
study;  however,  in  the  systematic  development  of  fragment  density  calculation, 
it  became  necessary  to  include  stack  models. 

1.2  Method  of  Approach 

A  review  of  the  existing  experimental  data  on  far-field  fragment 
densities  summarized  the  body  of  data  in  existence.  The  existing  techniques 
for  calculating  the  far  field  fragments  were  reviewed.  Using  these  techniques 
as  a  springboard,  an  approximate  technique  for  evaluation  of  the  far-field 
fragment  densities,  assuming  a  known  source  on  a  unit  hemisphere  at  ground 
zero,  was  devised  to  rapidly  estimate  the  effects  of  the  models  that  were 
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derived.  Comparisons  with  the  more  sophisticated  calculations  for  a  single 
munition  using  measured  arena  data  as  the  init:-!  data  were  favorable. 

Models  for  tht  effect  of  an  igloo  on  the  initial  data  were  prepared. 

The  models  were  (1 l  that  the  fragment  from  the  munition  must  perforate  the 

earth  cover,  (2)  tf  at  the  fragments  pass  through  an  earth  cover  that  is  in  the 
process  of  breaking  up  and  that  a  certain  fraction  are  stopped  while  the 
remainder  pass  through,  and  (3)  that  the  crown  section  of  the  earth  cover  is 
blown  off  and  does  not  effectively  stop  any  of  the  upper  register  fragments. 

A  model  for  the  dependence  of  initial  hemisphere  data  for  a  stack  of 
munitions  was  prepared  by  considering  that  the  fragments  from  each  munition 
must  perforate  the  expanding  metal  shell  from  any  adjacent  munition  and 
further  must  penetrate  the  dense  ietonation  products  from  the  explosive 
filler  of  each  adjacent  munition. 

Comparisons  of  each  theory  were  made  with  experimental  data  and 
the  models  that  best  fit  the  data  were  selected. 

In  addition,  it  is  known  tnat  initial  fragment  velocities  and  mean 
distributions  of  fragments  from  accidental  detonations  differ  from  the 
arena  data.  A  parameter  study  was  conducted  to  evaluate  these  ef feces.  To 
better  fit  the  far  field  data,  consideration  of  another  parameter,  namely, 
the  effective  ballistic  density,  was  included  in  this  parameter  study. 

A  tentative  model  for  quantity-distance  was  developed  from  the  results 
of  this  study  and  typical  quantity-distance  curves  were  calculated  for  four 
munitions. 

1.3  Main  Results 

The  model  which  most  adequately  describes  the  effect  of  the  igloo  is 
that  the  crown  of  the  igloo  Is  blown  off.  Only  fragments  with  Initial 
elevation  angles  less  than  10°  are  retarded  (these  are  completely  stopped). 
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About  a  10%  reduction  in  fragment  densities  in  the  far-field  is  effected  by 
the  presence  of  an  igloo.  The  model  for  the  stack  shows  that  only  the 
munitions  in  the  surface  layers  contribute  significantly  to  the  far-field 
fragment  density. 

The  initial  velocity  for  the  fragments  should  be  taken  as  the  highest 
value  from  the  arena  data.  The  fragments  that  are  thrown  to  the  far-field 
must  have  a  higher  effective  ballistic  density  than  the  mean  value  given 
in  the  arena  data;  i.e.,  the  fragments  that  reach  far  distances  tend  to  be 
elongated  with  the  long  axis  aligned  with  the  trajectory.  Since  these 
projectiles  have  a  high  mass/presented-area  ratio,  a  more  conservative 
injury  criterion  in  terms  of  impact  than  any  of  those  commonly  applied  thus 
far  may  be  appropriate. 

Because  there  is  very  little  data  on  far-field  fragment  measurements 
and  because  the  experiments  were  not  systematically  planned  to  measure  the 
dependence  of  the  far-field  fragments  on  the  several  parameters  of  the 
typical  munition  storage  configuration,  and  extensive  experimental  program 
to  determine  this  is  highly  recommended.  Until  experimental  verification 
of  the  predicted  parameter  dependence  is  experimentally  verified,  the 
analytical  predictions  given  in  this  report  must  remain  hypothetical. 
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Chapter  II 

EXPERIMENTAL  AND  THEORETICAL  BACKGROUND 


Several  studies  on  far-field  fragment  distributions  from  the  detonation 
of  munitions,  both  theoretical  and  experimental,  have  been  conducted.  A 
summary  of  these  5-tudies  is  given  in  the  following  paragraphs. 

2.1  Experimental  Studies 

Experimental  determination  of  fragment  distribution  in  the  far-field 
has  been  accomplished  for  a  limited  number  cf  munition  stores.  Because  of 
the  costs  and  time  for  detailed  experimental  determination  of  the  far-field 
fragment  hazard,  these  experiments  are  few  in  number  and,  since  in  most 
inscances  the  primary  purpose  of  the  experiment  was  blast  effects,  evalu¬ 
ation  of  structural  damage, or  detonation  communication,  the  fragment  study 
was  a  secondary  goal. 

In  1945,  the  ANESB  undertook  a  series  of  detonations  of  various  munitions 
contained  within  standard  Arn\y  and  Navy  storage  igloos  (ANESB  1947).  A 
series  of  eight  detonations  of  large  stacks  was  conducted;  blast  wave 
propagation  and  structural  damage  to  barracks  by  blast  were  the  primary 
objective.  Detailed  fragment  measurements  were  made  on  only  one  shot 
containing  425  MK33, 1100-lb  bombs  in  a  Navy  igloo.  There  were  250,000  lbs. 
of  50/50  Amatol  in  the  stack.  Figure  2-1  shows  the  fragment  density 
(in  number  of  fragments  per  square  foot)  vs.  range  along  the  line  where  the 
maximum  fragment  densities  occurred.  Qualitative  descriptions  of  fragments 
from  the  remaining  experiments  and  from  the  subsequent  model  igloo  tests 
indicated  that  this  fragment  pattern  was  replicated. 
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A  series  of  four  stack  detonations  in  revetments  was  carried  out  in 
the  Big  Papa  tests  (Peterson  et.  al .  1968).  The  first  two  detonations  were 
rectangular  stacks  of  Ml  17,  750-pound  bombs  and  M66A2  2,000-pound  bombs,  each 
bomb  being  prmed.  The  third  shot  was  a  similar  mixed  stack,  but  in  a 

hexagonal  stack  arrangement.  Fragment  densities  were  measured  on  those 

2 

shots;  in  the  first  two  the  collection  areas  were  90,000  ft  areas,  with 
only  three  collection  areas  within  3500  feet  from  the  stack.  On  the  third 
shot  the  collections  were  smaller  and  many  more  were  placed  from  1000  feet 
to  2500  feet.  The  fourth  shot  was  a  stack  of  Ml  17  750-pound  bombs  in  a 
rectangular  array  with  one  bomb  primed.  The  rest  detonated  by  blast  or 
fragment  effects.  No  far-field  fragment  distributions  were  measured  in  this 
shot.  Figure  2-2  shows  the  fragment  density  for  the  first  two  shots  along 
a  line  normal  to  the  axes  of  the  bombs. 

The  remaining  experimental  data  on  fragments  comes  from  close-in 
measurements  of  fragments.  Several  military  engineering  installations  have 
measured  fragment  properties  in  arena  tests  in  which  one  horizontal  munition, 
primed,  is  detonated  with  the  fragments  intercepted  a  short  distance  away. 

As  a  function  of  azimutn  angles  (measured  from  the  nose),  an  average  initial 
fragment  velocity  and  fragment  distributions  with  mass  are  measured.  The 
arena  results  are  summarized  in  the  JMEM  manuals  (JMEM  1970). 

Draper  and  Watson  (1970)  report  some  modified  arena  data  for  stacks 
of  munitions,  in  particular  a  sequence  of  experiments  on  stacks  of  155mm 
TNT-filled  projectiles.  Their  interim  report  on  raw  fragment  collection 
data  gives  a  qualitative  indication  of  how  fragment  mass  distributions  and 
initial  fragment  velocities  change  as  the  stack  size  changes.  They 
tested  stacks  3x3x1,  5x4x1,  6x6x1,  12x6x1,  and  40x6x1. 


GENERAL  AMERICAN  RESEARCH  DIVISION 


6 


FRAGMENT  DENSITY 


Far-field  fragment  distributions  were  measured  from  small  stacks  of 
Ml  17,  750-pound  tri tonal -filled  bombs  in  a  test  conducted  at  NWC  (Feinstein 
and  Nagaoka  1970).  The  stacks  were  2x3x1  and  5x3x1  (length  x  height  x  depth). 
All  bombs  were  primed.  Fragment  and  mass  densities  were  measured  from  500 
feet  to  1500  feet  from  the  stacks  on  the  nose,  side  and  base  directions. 

Figure  2-3  shows  the  fragment  density  in  the  side  direction  from  these  two 
shots. 

A  10x100x1  stack  of  155  mm,  TNT  loaded,  projectiles  was  detonated  at 
Yuma  (Feinstein  and  Nagaoka  1970)  with  fragment  and  mass  densities  measured 
in  the  nose,  base  and  side  directions.  One  projectile  was  primed. 

A  large  stack  (13,696)  of  155  mm  projectiles  contained  inside  a  standard 
Army  igloo  was  detonated  in  Eskimo  I  (Weals  1973).  The  projectiles  had  TNT 
as  a  filler.  Their  axes  were  vertical  and  approximately  20  of  these  projec¬ 
tiles  were  primed.  Figure  2.4  shows  the  measured  fragment  densities  along 
several  lines  from  this  experiment.  Figure  2-5  shows  the  plan  view  of  the 
experiment  and  defines  the  rays  relative  to  the  igloo  structure. 

2.2  Theoretical  Studies 

Very  few  theoretical  predictions  of  far-field  fragment  patterns  from 
munitions  have  been  made.  Zaker,  et  al.,  (1970)  generated  a  method  for 
calculating  far-field  fragment  densities  from  the  arena  data  for  a  single 
munition.  For  a  given  interval  of  fragment  mass  and  an  initial  velocity, 
(given  as  a  function  of  azimuth  angle  measured  from  the  bomb  nose),  they 
compute  the  ranges  wherein  the  fragments  will  fall  (for  all  possible 
elevation  angles).  In  each  increment  of  range,  the  number  of  fragments  over 
all  mass  intervals  is  summed  and  the  fragment  densities  are  then  calculated. 
They  present  predictions  for  far-field  fragment  densities  for  stored  single 
munitions.  Feinstein  (1972)  extended  this  calculation  with  an  interpolated 
mass  distribution  for  the  larger  mass  intervals.  Schreyer  and  Romesberg 
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(1970)  (also  Romesberg  1971)  prepared  a  similar  program  to  calculate 
fragment  distributions  that  would  be  expected  in  the  Big  Papa  experiment. 

In  all  these  programs,  the  initial  data  on  the  fragments  was  taken 
from  the  arena  data.  In  Zaker's  and  Feinstein's  work,  the  distribution  over 
mass  was  based  in  tabular  form  from  the  arena  data,  while  Romesberg  and 
Schreyer  used  the  mean  fragment  mass. 

Feinstein  attempted  to  extend  the  single  munition  calculation  procedure 
using  an  angle  dependent  factor  for  his  initial  data  obtained  from  experimental 
arena  data  on  model  stack  munitions  (Feinstein  and  Nagaoka  1970).  Schreyer 
and  Romesberg  determined  an  effective  number  of  single  munitions  by  comparing 
the  single  munition  prediction  with  an  experiment  from  the  Big  Papa  series. 

2.3  Some  Comments 

To  summarize,  currently  available  information  consists  of: 

1)  theoretical  calculations  for  far-field  fragment  densities  from 
single,  unbarricaded  munitions,  and 

2)  several  experimentally-measured  far-field  fragment  densities  from 
stacked  munitions  in  a  variety  of  con  figurations. 

Parametric  representations  of  far-field  densities  from  stacked  munitions, 
with  or  without  an  earth  cover,  are  virtually  non  existent. 

Very  little  work  has  been  done  In  quantitatively  evaluating  the  far- 
field  fragment  hazards  from  accidental  detonations  of  stored  munitions.  The 
experimental  studies  are  far  too  sparse  to  allow  any  interpolation  or 
extrapolation  for  ranges  at  which  specific  fragment  densities  will  occur. 

The  various  experiments  vary  from  one  type  of  munition  to  another;  the 
orientation,  mode  of  ignition  and  protective  covering  vary  without  any 
systematic  pattern.  Only  In  the  fragmentary  (no  pun  intended)  results  of 
Draper  and  Watson  (1970)  is  there  any  systematic  investigation  of  the  effect 
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tmtLn 


«•**»**.-*>  - . 


of  varying  a  single  parameter  (in  this  case,  the  number  of  munitions  in  a 
stack  n  x  m  x  1).  Any  predictive  technique  requires  some  knowledge  of  how 
the  basic  physical  parameters  of  interest  (e.g. , initial  fragment  velocity) 
vary  with  the  parameters  that  describe  the  storage  configuration  (e.g. .number). 
Any  theoretical  calculation  of  far-field  fragment  densities  must  incorporate 
the  dependence  on  munition  storage  parameters  and  the  check  on  any  theory 
is  its  comparison  with  experimental  results.  With  the  current  data  avail¬ 
able,  these  comparisons  are  extremely  restricted. 
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Chapter  III 

BASIC  MATHEMATICAL  ANALYSES 


3. 1  Approximation  to  Solutions  of  the  Ballistic  Equations 

The  solution  of  the  ballistic  equations  to  obtain  the  range  of  a 
fragment,  given  its  initial  velocity,  elevation  angle, mass,  presented  area 
and  drag,  is  necessary  to  the  evaluation  of  the  fragment  densities  that  may 
be  expected.  The  differential  equations  that  describe  the  motion  and  tra¬ 
jectory  are  nonlinear  and  the  solutions  to  these  equations  are  either 
numerical  or  approximate.  The  equations  of  motion  are  solved  numerically 
and  then  approximations  to  these  solutions  which  facilitate  the  evaluation 
of  the  far- field  fragment  densities  are  devised.  Consider  the  motion  of  a 
fragment,  initially  propelled  with  velocity,  V0,  at  an  angle,  otQ,  with 
the  horizontal  plane.  This  motion  is  confined  to  two  dimensions,  i.e.,  no 
yaw  effects  and  no  cross  wind  effects  are  considered.  The  forces  acting 
on  the  fragment  are  drag  and  gravity.  Figure  3-1  depicts  the  trajectory 
and  defines  the  coordinates.  The  equations  of  motion  are 


dV 

x 

m  dt 
dV 

y 

m  dt 


V  V 
x 


mg 


where 


dX  .  v 
dt  ¥x 


dY 

Ht 


V  2  +  V  2 
x  y 


m  is  the  mass  of  the  fragment 
CQ  is  the  drag  coefficient 


(3-1) 
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p  is  che  air  density 

A  is  the  cross  section  area  of  the  fragment  normal  to  the 
trajectory 

g  is  the  gravitational  acceleration 

V  ,  V  are  the  horizontal  and  vertical  components  of  the 
x  y 

fragment  velocity 
t  is  the  time 

The  initial  conditions  are,  at  t  =  0 
X  =  Y  =  0 

V  =  V  COS  a 

X  0  0 

Vy  «  VQ  sin  aQ 

where  VQ  is  the  initial  fragment  velocity  and  aQ  is  the  initial  elevation 
angle. 


The  presented  area,  A,  can  be  represented  approximately  as  a  function 


of  the  fragment  mass  through 


A  =  (1) 

tc ' 


(3-2) 


where  k  is  the  ballistic  density.  For  a  tumbling  fragment,  this  area  A  is 
taken  as  the  mean  presented  area, 
where  k  is  the  ballistic  density. 

Define  a  parameter,  6,  through 


(3-3) 


Then  the  first  two  equations  of  (3-1)  can  be  written. 


avx  2/3 

■  3T  -  -  V 


m  ■  -  err'  Vy V  -  mg 


(3-1 1 ) 
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The  drag  coefficient  is  a  function  of  the  velocity  as  is  shown  in 
Figure  3-2.  For  fragment  velocities  exceeding  the  sound  velocity  in  air, 
the  drag  coefficient  is  essentially  constant  (C^  ^  1.28),  while  for  velo¬ 
cities  below  the  sound  velocity  in  air,  it  has  another  constant  value, 

(Cg  %  1.08).  Only  when  the  fragment  velocity  is  close  to  the  sound 
velocity  does  the  drag  coefficient  vary  significantly.  The  fragments  under 
consideration  have  an  initial  velocity  several  times  the  sound  speed  in  air; 
most  of  the  trajectories  of  these  fragments  are  accomplished  while  the  fragment 
velocities  are  supersonic.  Thus,  for  simplicity  of  analysis  and  scaling, 
assume  the  drag  coefficient  to  be  constant  with  its  supersonic  value, 

CD  =  1.28. 

The  equations  of  motion  are  to  be  solved  from  t  =  0,  up  to  a  time  when 
y  returns  to  zero,i.e.,  when  the  fragment  trajectory  intersects  a  horizontal 
plane  through  the  source.  The  value  of  x  at  that  time  is  the  range,  R,  of  the 
fragment.  The  velocity,  V^,  is  the  impact  or  terminal  velocity  and  the 
trajectory  makes  an  angle,  a^,  with  the  horizontal. 

What  is  desired  are  the  solutions 

R=  fl(V0,  V  m) 

V  V  V  V 

af=  W  “o’ 

Introduce  the  new  independent  variable,  S,  the  arc  length  along  the 


(3-4) 


(3-5) 
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ORAG  COEFFICIENT 


Figur*  3-2.  DRAG  COEFFICIENT  FOR  A  FRAGMENT  AS  A  FUNCTION  OF  FRAGMENT 

VELOCITY  . 


dx  _  !x 

dS  '  V 

dY  _  ^ 

d5'  v 

The  initial  conditions  are:  at  S  =  0 
V  =  V  COS  a 


1/  sin  a 
o  o 


y 

x  =  y  =  o 

The  equations  are  put  into  a  dimensionless  form  through  the  introduction 


of  a  length  scale,  L,  and  velocity  scale,  V£ 

1/3 

L  =  - 
B 

The  equations  of  motion  can  be  written 

.  .  V 

ds  vx 


(3-6) 


<K 

di1  =  "  Vy 


£ 

V 


dx 

d? 


v 

_  X 


&  =  3L 

ds  v 


(3-7) 


.2 


2  ,  2 
=  v  +  v 
X  y 

t!«l 

BV  2 


X  Y  %  ¥v  J  S 

with  x  =  p  y  *  ,  vx  »  y-,  vy  «  ^  ,  and  s  =  j- 

The  parameter,  e,  is  the  ratio  of  the  terminal  free  fall  velocity  of 
the  fragment  to  its  initial  velocity. 
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The  initial  conditions  are,  at  s  =  0 


Vx  =  cos  a0 


Vy  =  Sin  ao 


x  =  .y  =  0 

In  this  form,  the  equations  are  to  be  solved  as  a  function  of  s  until 
s=s^  where  y(s^)  =  0.  The  values  of  x  and  v  at  this  value  of  s  define 
the  range  r,  and  terminal  velocity. 


r  =  x^)  =  f^.e) 
vVvx4(sl)  +  vy2(sl>  =  f2(V£> 


a^=  tan 


-1 


Vfi> 


W 


=  V“o,e) 


The  three  functional  forms  f^ ,  fg,  and  f^,  are  determined  by  numerical 
integration,  using  a  Runge-Kutta  technique. 

The  parameter,  e,  is  a  function  of  the  initial  fragment  velocity,  its 
mass  and  its  ballistic  density.  From  the  arena  data  for  these  fragment  para¬ 
meters. the  range  of  e  which  is  applicable  to  this  problem  is  determined.  The 
fragment  masses  range  from  .02  lb  to  .3  lb  and  the  initial  velocities  from 
2000  ft/sec  to  10,000  ft/sec.  The  ballistic  density  datu  is  given  by  a 


single  value,  determined  from  the  mean  value  of  the  average  presented  area. 


This  mean  value  is  typically  600  grains/inch  . 


Using  this  value  of  <,  and  letting  p  =  525  grains/ft  . 

J/6 


e  =  137.04 


m 


% 


3  =  1.7146  x  10 


-3 


Figure  3-3  shows  curves  of  constant  e  as  a  function  of  VQ  and  m. 


The  parameter  e  varies  between  .01  and  .04  for  the  dominant  portion 
of  the  arena  data.  The  curve  with  e  =  0.02  goes  through  the  middle  of  the 
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range  of  the  arena  data.  We  will  be  concerned  with  slightly  larger  frag¬ 
ments,  (e.g.,  1.0  lb  at  initial  velocities  near  5000  ft/sec),  for  which 
e  =  0.02  is  representative  of  these  larger  fragments  also. 

If  we  take  k  twice  its  mean  value,  then 
e  as  137.04  21/3 

Vo 

and  the  range  of  e  shifts  to  from  .01  to  .06. 

The  non-dimensional  range  and  terminal  velocities  were  calculated 
for  values  of  e  from  .005  to  0.10.  The  results  are  shown  in  the  next 
three  figures.  Figure  3-4  shows  r  as  a  function  of  aQ  for  various  values 
of  e.  Figure  3-5  shows  versus  aQ  and  Figure  3-6  shows  versus  aQ. 

The  range  versus  elevation  angle  curves  show  very  little  variation  in 
shape  in  the  range  0.01  <  e  <0.06.  The  maximum  of  the  curve  shifts  slightly 
from  18°  at  e  =  0.01  to  22.5°  at  e  =  0.6.  Thus,  if  these  curves  are  normal¬ 
ized  such  that  their  maxima  are  unity,  we  find  the  curves  almost  identical. 
From  Figure  3-3,  note  that  the  curve  for  e  =  .02  goes  through  the  middle 
of  the  arena  data.  Deviations  in  shape  of  the  range  curves  from  this 
are  minimal.  Thus,  the  ballistic  curves  for  e  =  .02  are  selected  to  describe 
all  ballistic  parameters. 

Note  the  changes  in  shape  of  the  and  curves  are  quite  small  in 
the  range  0.1  <_  e  <_  .06. 

Figure  3-7  shows  the  maximum  of  the  range  versus  elevation  angle  curve 
as  a  function  of  e.  This  curve  can  be  approximated  to  a  high  degree  of 
accuracy  by 

rmax  *  0,5140  "  1,0358  1n  e  (3-9) 

We  use  then  an  approximation  of  the  form 

r<Vc>  =  rmax  W  f<“o>  (3-,0) 
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Fiflur#  3-4,  BALLISTIC  TRAJECTORY, 

DIMENSIONLESS  RANGE  VS  ELEVATION  ANGLE 
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Figur.  3-5,  IMPACT  ANGLE  VS  ELEVATION  ANGLE. 


IMPACT  VELOCITY 


2.0 


VT  -  TERMINAL  FREE  FALL  VELOCITY . 
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Figure  3-6,  IMPACT  VELOCITY  VS  ELEVATION  ANGLE. 
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MAXIMUM  NON  DIMENSIONAL  RANGE 


Figurt  3-7,  MAXIMUM  NON  DIMENSIONAL  RANGE  VERSUS  €  -  RATIO 
IMPACT  VELOCITY /TERMINAL  VELOCITY. 
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where  f(aQ)  is  determined  from  the  curve  for  e  =  0.02.  After  some  trial 
and  error,  an  approximate  fit  for  f(aQ)  was  obtained. 

f(«0)  =  ^  (1-5)[1+A(5  -  y)2]  =  h(5)  (3-11) 

with 

A  =  3.25,5  <  j 
A  =  2.00,5  >  l 
*  =  sin  aQ 

Figure  3-8  shows  the  approximation  to  f(aQ)  versus  elevation  angle. 

The  fit  is  within  .1%  for  7.5°  <  a  <40°. 

o  — 

Other  approximations  are 

of  ~90  [1  -  exp  (-  -|^]  (3-12) 

~TS  (“  -  1  -9)  exp  (-  -|y)  +  1.6  exp  (-  -^-) 

-  0.7  exp  (-  -y^-)  +  1 
3.2  Calculation  of  the  Fragment  Density 

An  approximate  technique  for  the  rapid  estimation  of  fragment  densities 
in  the  far-field  as  a  result  of  a  detonation  of  a  munition  store  is  derived 
in  this  section.  This  approximate  procedure  utilizes  the  approximate  forms 
for  the  ballistic  parameters  developed  in  the  previous  section. 

Consider  ranges  which  are  sufficiently  large  so  that  the  munition 
store  can  be  treated  as  a  point.  Figure  3-9  shows  the  coordinate  system 
that  we  use.  The  point  where  we  will  compute  the  fragment  density  is  defined 
by  a  range,  R,  measured  from  the  center  of  the  munition  store  and  an  azimuth 
angle,  +,  measured  from  some  reference  axis  on  the  ground  surface.  When 
the  munitions  detonate,  the  fragments  at  the  source  are  described  as  follows. 
On  a  unit  hemisphere  surrounding  the  origin,  the  number  of  fragments  per  unit 
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Figure  3-8,  APPROXIMATION  TO  NONDIMENSIONAL  RANGES  VS 
ELEVATION  ANGLE,  FIT  TO  €  •  0.02. 
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V  ( a0<4>  ) 


Figurt  3-9,  ILLUSTRATION  OF  COORDINATE  SYSTEM  AND  UNIT 

HEMISPHERE  USED  FOR  CALCULATION  OF  FAR -FIELD  FRAGMENT  DENSITIES- 
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am 


solid  angle  in  a  given  mass  increment  are  known  as  a  function  of  azimuth 
angle  and  elevation  angle.  Also  the  initial  velocity  of  the  fragments  is 
known  as  a  function  of  these  angles.  Thus  at  fixed  aQ  and  <j>,  there  is 
a  probability  density  function  f(m,  aQ,  $)  such  that 
f (m,  aQ,  $)  dm 

represents  the  number  of  fragments  per  unit  solid  angle  with  mass  between 
m  and  m+dm.  Further  define  a  cumulative  distribution  function 


F(m,  a0,<j>)  =  f°  f(m,  aQ,  4)  dm  (3-14) 

m 

This  defines  the  number  of  fragments  per  unit  solid  angle  with  mass 
greater  than  m. 

In  an  infinitesimal  element  of  azimuth,  A$,  the  number  of  fragments 
that  fall  beyond  R  is  calculated  by  integrating  over  a  slice  of  the  unit 
hemisphere.  Denote  this  number  by  N (R ,«j>) . 


N(M) 


/  COS  aQ  daQ 


( R  *  aQ» 


f(m)dm 


A4> 


(3-15) 


The  lower  limit,  m^R,  ao>4>)  is  the  mass  of  the  fragment  that  just  reaches 
R  given  the  elevation  angle  aQ  and  the  initial  fragment  velocity  at  aQ  and 
4.  All  fragments  (with  the  same  $  and  aQ)  which  have  masses  larger  than  m^ 
(R,  aQ,  4)  will  fall  at  larger  ranges  than  R.  (Figure  3-10  illustrates  the 
domain  of  integration  in  the  aQ,  m  plane.)  By  considering  N(R,  4)  and 
N ( R  +  aR,  4)  for  an  infinitesimal  aR,  the  number  of  fragments  in  the  small 
area  RaRa$  is 


n(R,  AR,  ♦)  =  ARA*  (|^)  (3-16) 

The  fragment  density  at  R  is  then 

q(M)  *  (3*17) 
♦ 

Using  the  expression  for  N  from  Equation  (3*15), 
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■ |  tt/2 

q(R»  +  )  =  I  /  C0S  a0  da0  f(ml(R*ao,4>))(3R)a0,4) 

The  partial  derivative  (|§)  can  be  evaluated  from  the  expression 


for  the  range 


1/3  1/6  n  1/2 

R  =  j  h(5)  {.5140  -  1.0358  In  (|)  }  (3-18) 


a0»4>»m-| 


1/6  "  1/2 

{-.0039  -  1.0358  In  [p-  (J)  ]} 

Vo  3 


(3-19) 


To  evaluate  this  integral,  an  approximation  for  F(m,a0,$)  is  required.  A 
reasonable  approximation  to  the  arena  data  is  the  Thomas  (1944)  approximation 

F(m »a0 »4>)  =  N0(a0»$)  exp  {-  (m/y(a0>(j>))1^3} 

where  N0(a0*4>)  and  u ( aQ , )  are  parameters  indicating,  respectively,  the 
total  number  of  fragments  per  unit  solid  angle  of  aQ  and  4>Qand  one-sixth 
the  mean  fragment  mass  in  the  direction  given  by  aQ  and  <j>. 

From  tnis  expression  for  F(m,a0>4>),  it  follows  that 

=  ~T/T~a£~2'/3  exp  <-(m/ii(a0,4»))1/3}  (3-20) 

3u  (a0  »4>)m  ' 

Then,  with  £  =  sin  aQ  as  the  integration  variable, 

1  dcN0(a0,t)  exp  {-[m|(R,a0,«|>)/p(a0,4)]1;/3} 

q(R’^  R  o  (u(a  ,0)}1/3h(€)  {-.0039  -  1.0358  ln(^-/6(J)1/2)} 

v0  p 


(3-21) 


If  the  source  is  isotropic,  the  expression  simplifies  somewhat. 


,I<R)  =  J'3  l 


1  dt  exp  {-(m^R.aJ/u  )l/3} 


h(c)  j^- . 00 39  -  1.0358  ln(^  (*)  )j 


(3-22) 
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This  approximate  calculation  is  checked  by  calculating  *he  far- field 

fragment  density  in  the  side  spray  direction  for  a  single  munition  with  its 

axis  horizontal.  The  values  for  y  and  V  were  obtained  from  the  arena  data 

o 

by  averaging  these  numbers  in  the  range  60°  <  ^  <  120°,  where  ip  is  measured 
from  the  axis  of  the  munition  extended,  ^  =  0°  denoting  the  nose  end  of  the 
munition.  Table  3-1  show?  the  values  used  in  this  calculation. 

Table  3-1 


INITIAL  FRAGMENT  PARAMETERS  -  SIDE  SPRAY 
DIRECTION  -  AVERAGED  FROM  ARENA  DATA 

w  / Fragments^  v  /ft  x 

Munition  Explosive  Msteradi an/  vo  ^sec^  yi(lbs) 


155  mm  shell 

TNT 

500 

4000 

.010 

660 

175  mm  shell 

Comp  B 

3000 

5000 

.00429 

660 

500  lb  bomb 

H-6 

4500 

7000 

.00286 

590 

750  lb  bomb 

Tri tonal 

6000 

8000 

.00429 

590 

In  the  side  spray  direction,  there  is  no  dependence  of  any  parameters 
on  either  aQ  or  $,  and  the  integral  takes  on  the  isotropic  form  (Eq.  3-22). 
The  values  of  q  for  these  four  munitions  were  calculated  by  numerically 
evaluating  this  integral.  The  results  are  shown  in  Figures  3-11  through 
3-14.  Feinstein's  (1972)  calculated  values  of  q  in  the  direction  ip  -  90°, 
obtained  by  numerical  integral  of  the  arena  data  over  the  unit  hemisphere 
is  shown  for  each  case.  The  two  models  have  the  same  general  trend  and 
agree  throughout  the  range  500  ft-2000  ft  to  within  a  factor  of  2. 

An  estimate  of  the  error  induced  in  the  calculation  of  q  when  the 
ballistic  range  approximations  of  the  above  are  used  was  obtained  by  com¬ 
paring  the  solution  for  the  500  lb  bomb  wifn  a  calculation  using  the 
tabular  values  of  range  from  the  numerical  solution  to  the  equations  of 
motion.  The  error  in  q  was  from  0.15%  to  0.50%.  The  error  was  primarily  at 
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FRAGMENT  DENSITY 


E. 


Figure  3-11,  SINGLE  MUNITION  FRAGMENT  DENSITY  VERSUS  RANGE  MI07 
135mm  PROJECTILE  (TNT  LOADED  SIDE  DIRECTION) 
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Figure  3-13,  SINGLE  MUNITION  FRAGMENT  DENSITY  VERSUS  RANGE 
Mil  7  750-  lb.  BOMB  (SIDE  DIRECTION  ). 
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FRAGMENT  DENSITY 


Figure  3*14,  SINGLE  MUNITION  FRAGMENT  DENSITY  VERSUS  RANGE 
MK  82  500-lfo.  BOMB  (SIDE  DIRECTION) 
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very  small  aQ(0.  <  aQ  <  5°)  or  large  aQ(a  >  50°),  where  the  integrand  is 
small . 

From  an  examination  of  the  integrand  for  the  approximate  integration 
over  aQ,  a  very  simple  and  crude  approximation  to  q  may  be  obtained. 

q~  jj  exp  (-BR)  (3"23) 

where  A,  and  B  are  constants  for  each  weapon  (Table  3-2). 

Table  3-2 

CONSTANTS  A,  B  FOR  APPROXIMATION  TO  q 
Munition  A  (No/ft)  B  (  x  103  ft  ^) 


155  mm 

projectile 

.8177 

1.893 

175  mm 

projectile 

6.103 

2.380 

500  lb 

bomb 

10.32 

2.711 

750  lb 

bomb 

11.68 

2.309 

This  crude  approximation  agrees  with  the  integral  to  within  10%  for 
500  <  R  <  2000  ft. 

This  chapter  has  described  a  rapid  technique  for  estimating  fragment 
densities  and  mass  densities  in  the  far-field.  Further,  the  results  of 
this  approximate  procedure  compare  favorably  with  Feinstein's  (1972)  computer 
code  predictions.  The  basic  parameters  that  describe  the  single  munition 
are  the  initial  fragment  velocity,  the  ballistic  density  of  the  fragments 
and  the  form  of  the  probability  density  function  giving  the  number  of 
fragments  in  any  mass  interval.  When  the  more  complicated  problem  of 
munition  stacks  and  stacks  within  igloos  is  attached,  the  obvious  place 
to  start  is  to  assess  the  changes  in  the  source  parameters  induced  by  the 
more  complicated  configuration. 
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Chapter  IV 

ANALYTICAL  MODELS  FOR  STACK  AND  COVER  EFFECTS 
ON  FAR- FIELD  FRAGMENT  DENSITIES 


4.1  Introduction 

In  this  chapter,  several  models  are  described  which  attempt  to  explain 
the  effects  of  the  cover  and  the  stack  configuration  on  the  far-field 
fragment  density.  There  are  several  common  features  in  these  models  which 
will  be  qualitatively  summarized  here;  more  detailed  exposition  of  the 
features  will  be  found  in  the  appropriate  sections. 

The  initial  fragment  distribution  and  other  initial  parameters  when  a 
single  munition  is  initiated  under  the  influence  of  a  proper  fuze  are  well 
defined  through  physical  measurements.  When  a  stack  of  the  munitions  is 
Initiated  at  a  single  munition  and  the  remaining  munitions  detonate  by 
sympathetic  detonation,  induced  by  either  blast  or  fragment  impact,  the 
initial  fragment  distributions  and  initial  velocities  from  any  single  muni¬ 
tion  in  the  stack  may  be  drastically  altered.  For  example. Feinstein  and 
Nagaoka  (1970),  Feinstein  (1972a,  1972b)  found  in  the  detonation  of  1000 
(10  x  100  x  1)  155  mm  TNT-loaded  projectiles,  many  large  fragments  none  of 
which  ever  appeared  in  the  data  from  any  arena  tests.  Similar  observations 
were  made  at  the  Eskimo  I  test  (Weals  1973).  Draper  and  Watson  (1970)  report 
altered  fragment  distribution  patterns  and  increased  initial  fragment  velo¬ 
cities  in  the  nose  and  base  directions  from  their  stack  experiments. 

The  second  feature  is  that  fragments  from  one  munition  may  interact  with 
fragments  from  other  munitions,  the  blast  waves  from  other  munition  or  a 
piece  of  the  cover  structure.  Its  velocity  will  change  and  its  trajectory 
may  be  altered. 
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Consider  the  simultaneous  detonation  of  adjacent  stacked  munitions. 
Fragments  emanating  from  any  one  weapon  will  collide  and  interact  with 
fragments  from  all  the  other  weapons  in  the  stack.  Each  collision  will 
result  in  a  change  in  the  speed  and  direction  of  each  fragment  involved  in 
the  collision.  Hence,  the  overall  effect  of  the  collision  processes  will 
be  to  alter  the  speed  and  direction  of  all  fragments  involved  in  the  colli¬ 
sions;  i.e.,  the  collision  process  will  alter  the  mass  and  velocity  distri¬ 
butions  of  the  fragments  emanating  from  any  one  weapon.  Furthermore,  the 
altered  distributions  of  weapon  fragments  will  depend  on  the  location  of  the 
weapon  within  the  stack  and  the  geometry  of  the  stack. 

Now  consider  another  related  problem,  that  of  a  weapon  detonated  with¬ 
in  an  earth-covered  igloo.  The  blast  wave  from  the  detonation  causes  the 
earth  cover  to  move  outwards.  Shortly  thereafter,  fragments  from  the 
weapon  interact  with  the  earth  cover.  The  situation  is  similar  to  that 
above  in  that  the  overall  effect  of  the  'interaction  is  also  to  change  the 
mass  and  velocity  distributions  of  the  fragments  from  the  weapon. 

It  should  be  noted  that  the  most  complex  problem  to  be  considered,  that 
of  a  stack  of  munitions  detonating  within  an  igloo,  can  be  reduced  to  the 
two  problems  discussed  above  by  adopting  the  following  approach.  Essentially 
all  of  the  fragment-fragment  interactions  occur  within  a  region  limited  to 
the  immediate  neighborhood  of  the  weapons.  Furthermore,  it  is  reasonable 
to  assume  that  fragments  from  any  one  weapon  can  significantly  interact  only 
with  fragments  from  weapons  in  proximity  to  each  other.  The  altered  frag¬ 
ment  mass  and  velocity  distributions  for  this  group  of  munitions  can  be 
transposed  to  an  origin.  In  effect,  a  group  of  stacked  weapons  can  be 
replaced  by  a  single  equivalent  weapon  that  interacts  with  the  earth  cover 
as  mentioned  earlier. 
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4.2  Fragment-Fragment  Interactions 

Since  the  largest  number  of  fragments  from  a  single  munition  are  ejected 
to  the  side,  consider  only  this  side-spray  and  idealize  the  munition  to 
infinite  cylinders.  When  a  single  weapon  detonates,  the  cylindrical  shell 
breaks  up  such  that  all  fragments  leave  at  the  same  initial  velocity.  This 
initial  velocity  and  the  mass  distribution  of  the  fragments  are  considered 
independent  of  angle  and  are  known.  Because  the  region  of  concern  is  the 
immediate  neighborhood  of  the  weapon,  the  effects  of  gravity  may  be 
neglected,  so  that  the  trajectories  are  straight  lines  given  by: 

S  =  j  In  (1  +  em1/3  VQt)  (4-1) 

It  is  seen  that,  for  a  given  time,  there  is  a  maximum  distance  beyond 
which  no  fragments  will  be  found,  and  that  more  massive  fragments  travel 
longer  distances. 

The  mass  distribution  of  fragments  ejected  at  any  angle  is  assumed 
given  as 

1  /3 

f(«)  -  j  exp  {-  (m/w)1/3}  (4-2) 

where  y  =  (average  mass)/6  =  known  constant 

Choose  a  distance,  S  =  d,  and  a  time  t.  Let  m  be  the  zero  of  eq.  (4-1). 
The  fraction  of  total  fragments  of  mass  greater  than  m: 

tt®-  =  _/  f(m)  dm  =  exp  {-  (m/y)1^3}  (4-3) 

"o  in 

where  N(rh)  =  number  of  fragments  of  mass  greater  than  in. 

Nq  -  total  number  of  fragments 

•  -*  -.-4  4  i  ♦  ♦ 
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I  1 1  I  r 


The  number  of  fragments  located  at  a  distance  greater  than  d  at  a 


given  time,  t: 


N(d;  t)  =  N(m)  =  NQ  exp  {-  (m/p)1/3}  (4-4) 

Another  expression  for  N(d;  t)  is  simply 

N(d;  t)  =  /  p*  dV  (4-5) 

V 


where  p*  =  number  of  fragments  of  mass  iii  per  unit  volume 
In  polar  coordinates 

S 

N(d;  t)  =  2tt  /  max  p*  S  dS  (4-6) 

d  ' 

Equating  eqs.  (4-4)  and  (4-6)  leads  to  p*  (d;  t,  iii),  the  density  of 
fragments  of  mass  iii,  at  time  t,  at  distance  d: 


p*(d;t ,m) 


m  u  / 


S  =  d 


exp  {-(m/p)1/3} 


(4-7) 


The  derivative  is  found  from  eq.  (4-1).  Eq . (4-7 )  yields  fragment  densities 
only  in  the  immediate  neighborhood  of  the  weapon. 

Now  consider  two  munitions  with  their  axes  parallel  which  are  deton¬ 
ated  at  the  same  time.  If  no  collisions  between  fragments  are  considered, 
fragment  densities  from  two  weapons  would  be  simply  the  superposition  of 
the  fragment  densities  from  each  weapon  individually.  However,  if  fragments 
from  both  weapons  occupy  the  same  region  of  space  at  the  same  time,  there 

will  be  a  finite  probability  of  collision.  We  calculate  the  fraction  of  a 
differential  volume  element  (at  a  given  distance  and  time)  occupied  by 
fragments  from  each  munition  and  set  the  probability  of  collision  propor¬ 
tional  to  their  product. 
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Refer  to  the  collision  geometry  of  Figure  4-1.  Then  the  probability  of 


collision,  P(d^,  d2;  t),  is  given  by 


P(dr  d2i  t)  = 


aV]  aV2 
(AV)2 


(4-8) 


AV^  =  volume  occupied  by  fragments  of  mass, 
m,  from  weapon  1 

a V2  =  volume  occupied  by  fragments  of  mass, 
m,  from  weapon  2 

aV  •  f  **  r 

Now  —1  =  No.-  of  f r*^nts_o.f -  mass _mj  1  volUme  of  fragment  of  mass  m- 
A V  unit  volume  I  i 

AV.  mi 

or  w  ■  p/ldji  t.  m,)  - 

where  pf  =  fragment  mass  density 


Hence,  in  this  case, 

P(dj ,  d2;  t)  =  ml )]  [p2*^2*  ^2^]  *^1  (4-10) 


Equation  (4-10)  yields  the  probability  of  collision  between  a  fragment  of 
mass  m1  and  a  fragment  of  mass  m2  occurring  at  the  point  (d-j,  d2)  at  time  t. 

Obviously,  the  fragments  which  collide  are  the  only  ones  which  undergo 
velocity  transfunmations ;  the  fragments  which  do  not  collide  proceed  on  their 
original  trajectories.  Hence,  the  number  of  collisions  occurring  at  a  given 
point  in  space  must  be  known,  and  can  be  calculated  as: 


Ncol 1 ^dl ,d2;  ml 


m2>  "  [ 


N-j  ( d-j » d 2 it , in i 


>][ 


N2 (d.|  ,d2  jt  ,m2 


)]  [P(drd2;t 


»] 


(4-11) 

where  NCQ^  *  number  of  collisions  occurring  at  (d.|,  d2)  at  time  t 
between  fragments  of  mass  and  m2 
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(d^  »dg ; t ,rr,^ )  =  number  of  fragments  of  mass  m^  from  weapon  i 
at  (d.|,  dg)  at  time  t 

Number  of  fragments  of  mass  m^  from 
Now  assume  N^(d-|,  d2;  t,  m^ )  =<  weapon  i  present  at  (d-| ,  d^)  at  time 

t  in  the  absence  of  weapon  j 


(4-12) 


Number  of  fragments  of  mass  m.  from 
— <weapon  i  which  have  collided  prior  to  > 
reaching  (d-j ,  d2) 

Inherent  in  this  assumption  is  the  limitation  that  each  fragment  can 
undergo  only  one  collision.  N^  (d^,  d?;  t,  m^ )  can  be  expressed  in  terms 
of  p^*(d^;  t,  m^)  and  P(d^,  d2;  t),  and  is  therefore  a  known  function. 

Then,  the  number  of  collisions  for  all  space  and  time  are  known  for 
all  values  of  m. 

The  number  of  fragments  from  each  weapon  undergoing  collisions  at  any 
point  in  space  at  any  time  is  equal  to  the  number  of  collisions  occurring  at 
that  point  at  that  time,  and  is  known.  The  collision  process  is  illustrated 
in  Figure  4-2.  It  should  be  noted  that  the  case  illustrated  is  for  the 
situation  in  which  each  fragment  maintains  its  identity;  i.e.,  fragments 
do  not  split  due  to  the  collision.  The  equations  describing  the  collision 
of  two  fragments  are 

ml  Vlk  +  ”2  V2t  *  "l  Vlt  +  m2  \ 

"l  V]  "  ffl2  V2  "  Vl  *  m2  V2 

n  n  n  n 


e(V,  +  V2  )  =  V,  +  m2  V2 
•n  n  •  «n  »  n 


V  -  V, 
't  t 


2t  =  V 
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The  subscripts  1  and  2  refer  to  fragments  1  and  2;  the  unprimed  values 
refer  to  conditions  before  impact  and  the  primed,  after  impact.  The  first 
two  equations  describe  the  conservation  of  momentum,  while  the  third  describe 
the  energy  loss;  e  is  the  coefficient  of  restitution;  e  =  1  for  perfectly 
elastic  impact  and  e  =  0  for  perfectly  inelastic  impact;  in  general, 

0  <_e  £  1 . 

For  the  number  of  fragments  which  collide,  ,  the  results  of  the 
collision  can  be  determined.  Note  that  the  coefficient  of  restitution,  e, 
is  a  parameter  which  must  be  supplied;  this  allows  the  results  of  a  collision 
to  be  determined  assuming  perfect  elasticity,  perfect  plasticity,  or  cases 
in  between.  At  the  high  velocities  expected  in  this  problem,  e  ^  0.  From 
the  manner  in  which  the  results  of  the  collision  process  vary  in  space,  the 
number  of  collisions  occurring  in  space,  and  the  number  of  fragments  which 
do  not  collide,  the  new  mass  and  velocity  distributions- as  angular  functions 
can  be  determined. 

An  approximation  is  now  devised  for  the  estimation  of  the  interaction 
effects  of  two  munitions.  The  motion  of  a  typical  fragment  mass,  mp 
ejected  from  munition  1  at  the  angle  a-|  is  calculated  under  the  following 
assumptions.  We  ignore  the  changes  in  trajectory  direction  and  approxi¬ 
mate  the  additional  resistance  to  the  test  fragment  as  drag  effected  by 
a  continuous  cloud  of  fragments  emanating  from  munition  2.  First  the 
number  density  at(  ,,  o^Jat  time  t  from  the  second  munition  is  calculated. 
This  is  given  by  Eq.  (4-7).  With  the  mass  distribution  function  described 


N  1/3  1/ri 

P*(d2;t,m)  =  - exp  [-(m/u) l/JJ 

c  6*d2ul/J  d2  -  V  t  exp  (-eml/Jd2) 


with  <  V  t 
2  o 


(4-13) 
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This  density  is  calculated  for  the  four  munitions  of  interest  in  this  report. 
Figure  4-3  shows  P*  ,  the  number  density  of  fragments  versus  d,  for  the 
second  munition  (in  this  example,  the  Mil 7  750-lb  bomb). 

The  density  curves  show  a  slow  rise  from  the  arrival  of  the  first 
(and  largest  fragments)  followed  by  a  rapid  rise  near  the  propagation  of 
the  fragments  with  mean  mass,  6y.  Effectively  the  fragments  propagate 
outwards  as  a  shell.  The  propagation  characteristics  of  this  shell  are 
those  associated  with  the  mean  fragment  mass.  At  this  shell  radius,  the 
mass  density  rises  rapidly  to  a  relatively  large  value.  Hot,  condensed 
combustion  products  from  the  detonation  of  the  explosive  filler  are  con¬ 
tained  within  this  shell.  At  these  early  times,  the  hot  gases  have  not 
had  time  to  expand  and  their  density  is  very  nearly  that  of  the  solid 
explosive. 

Now  if  the  test  frc  is  ejected  at  such  an  angle  that  it  misses 
the  shell  of  fragments  from  the  second  munition,  the  probability  of  it 
colliding  with  any  fragments  from  the  shell  is  very  small.  If,  on  the 
other  hand,  the  fragment  intersects  this  shell,  the  probability  of  inter¬ 
action  is  almost  unity.  As  a  matter  of  fact,  many  such  collisions  may 
occur.  The  best  model  for  the  motion  of  the  fragment  is  that  the  fragment 
must  perforate  a  continuous  steel  shell  of  approximately  the  same  thick¬ 


ness  as  the  original  munition.  The  drag  force  for  this  is 

c  _  a  1/3  .,2  pcloud 
Fn  =  0m  V  - 

U  p  . 

ai  r 


(4-14) 


We  calculate  the  relative  reduction  in  the  fragment  velocity  as  it 


perforates  the  shell 


V —  =  exp  {-  (p  -  P  .Jh}  (4  15) 

ai  r 

Calculations  of  this  reduction  i velocity  for  fragment  mass  up  to 
0.01  lb,  and  original  munition  wall  thickness  up  to  .01  ft  shows  that  the 
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Figur*  4-3,  NUMBER  DENSITY  OF  FRAGMENTS  FROM  A  MII7  750-lb  BOMB  ALONG  A  RAY. 


velocity  reductions  are  greater  than  50%.  For  1  lb  fragments,  the 
velocity  reduction  is  about  10%. 

After  penetration  of  the  shell,  the  test  fragment  must  move  through 

the  hot,  condensed  combustion  products.  This  material  is  about  one  fourth 

as  dense  as  steel  and  the  length  of  traverse  of  a  fragment  within  the  hot, 

dense  gas  is  at  least  an  order  of  magnitude  larger  than  the  traverse  through 

the  steel  shell.  Thus  a  1-lb  fragment's  velocity  will  be  reduced  by  roughly 

40%  and  any  fragment  smaller  than  .25  lb  will  be  effectively  stopped. 

Thus,  the  test  fragment  from  one  munition  is  effectively  stopped  if 

it  intersects  the  shell  with  mass  u  expanding  from  the  second  munition.  The 

limiting  value  of  ,  for  which  fragments  of  mass  m,.  just  intercept  this 

shell  from  the  second  munition  can  be  readily  expressed. 

d? 

tan  a1  =  (4-16) 

where  (referring  to  Figure  4-1) 

d,  =  r  +  6m1/3  ln(l  +  6m1/3VQt) 
d2  =  r  +  £p1/3  1  n  ( 1  +  Svi1/3V0t) 
with  r  the  radius  of  the  munition. 

Bounds  on  a-j  are 

sin"^  (^)  <_  a-j  <  ^  when  m  >_  y 
where  a  =  inte munition  spacing 

For  test  fragments  with  m  ^  u  expand  the  logarithms  in  a  power  series 

r  +  V  t  -  Su1/3V  2t2 

tan  a,  = - - - T/tV?  (4"17) 

1  r  +  V  t  -  em  /Jvn  r 
o  o 

If  m  >_  u,  but  nearly  equal ,  ^  |- 

When  fragments  that  are  to  be  hurled  to  large  distances  are  considered, 
only  fragments  with  m  >  u  need  to  be  considered,  but  the  majority 
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of  those  are  not  too  much  larger  than  y.  Thus  the  value  of  j 
as  the  limiting  value  for  interception  of  the  fragments  of  interest  is 
a  good  approximation.  It  should  be  noted  that  the  assumption  a  -v  d-j  is 
required  for  this  approximation  to  be  valid.  If  the  intermunition  spacing  is 
significantly  larger,  the  limiting  value  of  d-j  will  be  correspondingly  smaller. 

The  conclusion  is  reached  t,sat  a  fragment  from  one  munition  will  be 
completely  stopped  by  the  fragments  and  the  gases  from  an  adjacent  munition 
if  the  (test)  fragment  intersects  any  fragment  with  mass  less  than  or 
equal  to  6y,  and  the  only  fragments  that  effectively  escape  from  the 
adjacent  munition  have  their  original  trajectories  outside  a  45°  angle 
from  the  line  joining  the  centers. 

The  case  of  three  or  more  adjacent  stacked  munitions,  while  more  com¬ 
plex  in  terms  of  mathemati cs ,  can  easily  be  detailed  on  a  conceptual  basis 
by  extending  the  approach  discussed  here.  Only  three,  or  perhaps  four, 
adjacent  weapons  need  to  be  considered,  since  fragments  from  a  weapon  will 
interact  significantly  only  with  fragments  from  weapons  adjacent  to  each 
other.  Stack  geometries  to  be  considered  for  3  and  4  weapons  are  shown  in 
Figure  4-4.  These  geometries  allow  determination  of  the  interactions  among 
fragments  from  the  weapons  shown  in  Figure  4-5. 

The  determination  of  the  mass  and  velocity  distributions  of  fragments 
from  a  multiple-weapon  stack  is  achieved  by  considering  portions  of  the 
stack  as  shown  in  Figure  4-5,  for  example,  replacing  them  by  equivalent 
single  weapons,  and  repeating  the  process  until  the  entire  stack  is 
replaced  by  an  equivalent  single  weapon. 

Consider,  for  example,  three  col  linear  munitions.  The  fragments 
that  emanate  from  the  center  munition  will  be  stopped  by  the  adjacent 
munition  on  either  side  in  an  angular  range  up  to  45°  from  the  line  joining 
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2  WEAPONS 


3  WEAPONS 


4  WEAPONS 


Figure  4-4,  STACK  GEOMETRIES  CONSIDERED  FOR  ANALYSI 


52 


Weapon  chosen  for  analysis 


Figure  4-5,  FRAGMENT-FRAGMENT  INTERACTIONS  DETERMINED  BY 

ANALYZING  THE  CONFIGURATIONS  OF  FIGURE  4-4  (EQUAL 

SPACING  BETWEEN  WEAPONS). 
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their  centers.  Each  end  munition  will  be  affected  as  in  the  case  of  just 
two  munitions.  It  is  a  very  simple  and  straightforward  process  to  build  up 
more  complicated  stack  models  based  on  this  approximate  model.  It  suffices 
to  note  that  no  fragments  will  escape  if  a  munition  is  completely  surrounded. 

All  fragments  from  munitions  within  the  stack  are  sufficiently  retarded 
by  their  neighbors  such  that  they  will  not  reach  the  far-field.  The  frag¬ 
ments  from  one  munition  in  the  surface  layer  from  45°  to  135°  will  be  thrown 
to  the  far-field.  Thus  the  following  model  for  the  effective  number  of 
munitions  in  a  stack  is  generated.  Let  N<-  be  the  number  cf  munitions  in  this 
side  layer  and  the  number  of  munitions  in  the  top  layer.  Then  the  far- 
field  fragment  density  is  approximately 

-,n  */2 

q(R)  =  Ns  r4  F  (a  ,R)daQ  +  /  FT(a  ,R)da0  (4-18) 

o  u/4 

where  F(aQ,R)  is  defined  by  the  approximation  defined  in  Chapter  III.  If 
both  layers  have  the  same  munition  orientation  relative  to  the  normals 
to  the  respective  surfaces,  we  can  approximate  the  integrals  as 

q(R)  =  (0.9  Ns  +  0.1  NJ)q1  (4-19) 

where  q^  is  the  far-field  fragment  density  from  a  single  munition. 

In  this  model  we  have  neglected  the  changes  in  direction  of  the  frag¬ 
ments  due  to  collision.  For  the  larger  fragments  that  would  interact  only 
with  the  larger  fragments  from  the  adjacent  munition,  the  probability  of 
collision  is  very  small.  One  other  feature  that  has  been  neglected  is 
that  the  blast  waves  from  two  adjacent  munitions  will  interact  and  the 
resulting  pressure  field  will  alter  the  trajectories  of  fragments.  The 
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blast  interaction  effects  are  minimal  in  the  direction  normal  to  the  line 
joining  the  munitions. 

4.3  Fragment-Cover  Interactions 

Look  next  at  the  effect  of  the  igloo  structure  on  the  fragments.  We 
assume  that  the  fragment  distribution  froir;  the  stacks  are  known  and  that 
they  can  be  represented  cn  a  unit  hemisphere  located  within  the  confining 
structure.  At  each  direction  on  this  unit  hemisphere,  the  initial  frag¬ 
ment  velocity,  the  mass  distribution  and  the  shape  of  the  fragments  are 
known.  Now  if  there  were  no  earth  covers,  it  would  be  simple  matter  to 
compute  the  far-field  fragment  density.  The  cover,  however,  will  retard 
or  stop  some  or  all  of  one  fragments.  In  me  next  few  paragraphs,  we  will 
construct  three  models  that  represent  possible  interactions  of  the  cover 
with  the  fragments. 

The  typical  earth-covered  igloo  structure  is  a  cylindrical  steel 
arch,  with  a  rectangular  floor  plan;  the  front  and  rear  wall  are  vertical. 

The  rear  wall  and  the  cylindrical  portions  are  covered  with  earth,  (Figure 
4-6  shows  the  cross  section  through  the  cylindrical  arch);  the  earth  cover 
has  a  minimum  of  two  feet  at  the  crown.  The  front  wall  is  of  reinforced 
concrete  and  has  doors  of  structural  steel  plate. 

When  the  stack  detonates,  a  fraction  of  the  energy  goes  into  breaking 
up  the  munition  cases,  a  fraction  into  kinetic  energy  and  the  fragments  to 
be  thrown  out,  and  a  fraction  remains  in  the  blast  wave.  This  last  fraction, 
which  remains  in  the  air  blast,  can  be  estimated  by  the  Fano  formula.  The 
blast  will  interact  with  cover  and  cause  it  to  deflect  and  break  up.  A 
certain  fraction  of  this  blast  energy  will  be  absorbed  in  the  breakup  of 
the  cover  and  in  the  kinetic  energy  of  the  cover  fragments,  (which  will 
consist  of  steel  arch  fragments,  earth  clods,  concrete  fragments  ana  pieces 
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CONCRETE  FLOOR 


CROSS  SECTION  THROUGH  ARMY  STANDARD  IGLOO  MAGAZINE 


of  the  door).  From  the  studies  of  the  blast  that  is  propagated  to  large 
distances  from  explosions  within  igloos  this  fraction  of  energy  absorbed 
by  the  cover  is  quite  small,  being  of  the  order  of  10%  of  the  blast  energy 
released  from  the  open  stack.  If  all  of  this  absorbed  energy  shows  up  as 
kinetic  energy  of  the  earth  cover,  the  typical  earth  cover  fragment  velo¬ 
city  will  be  of  the  order  of  1000-2000  ft/sec. 

The  assumption  is  made  that  in  the  case  of  a  stack  of  weapons  contained 
within  an  earth  cover  all  interactions  between  fragments  from  adjacent 
weapons  occur  substantially  before  any  fragments  interact  with  the  earth 
cover.  With  this  assumption,  the  stack  can  be  replaced  by  a  single  equivalent 
weapon  with  known  mass  and  velocity  distributions. 

Consider  a  cylindrical  shell  contained  within  a  semi -cylindrical 
earth  cover.  When  the  munition  is  detonated,  the  resulting  blast  wave  creates 
an  overpressure  on  the  inside  of  the  earth  cover,  causing  it  to  move  radially 
outwards  in  a  known  manner*,  see  Figure  4-7.  The  earth  cover  is  considered 
to  have  no  strength  in  the  tangential  direction.  As  the  earth  cover  moves 
out,  it  breaks  up  into  small  fragments  which  are  contained  between  an  inner 
and  an  outer  radius.  The  cover  fragments'  mean  motion  can  be  determined 
from  Newtons'  Law.  If  gravity  is  again  neglected,  all  interactions  between 
fragments  and  the  earth  cover  occur  along  radii;  i.e.,  the  only  result  of 
a  fragment-cover  interaction  is  to  alter  the  magnitude  of  a  fragment's 
velocity,  but  not  its  direction. 


*  Depending  on  the  relative  sizes  of  the  stack  and  the  igloo,  either  the 
blast  wave  or  the  fragments  will  arrive  at  the  earth  cover  first.  The 
effect  of  the  blast  wave  in  either  case  is.  considered  independent  of 
the  effect  of  the  fragments;  the  blast  wave  merely  specifies  initial 
conditions  for  the  fragment-cover  interactions.  The  following  analysis 
applies  in  either  case. 
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gure  4-7,  MOOEL  FOR  FRAGMENT- COVER  INTERACTIONS  (ASSUMING 
BLAST  WAVE  ARRIVES  AT  COVER  BEFORE  FRAGMENTATION 
WAVE  ). 
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Define  three  regimes:  pre-interaction,  interaction,  and  post-inter- 
action.  Only  air  drag  affects  the  fragments  during  the  pre-  and  post¬ 
interaction  regimes.  If  the  cover  is  considered  as  a  continuum  of  infinitesi¬ 
mally  small  particles,  the  retardation  force  between  the  fragment  and  the 
earth  cover  during  the  interaction  regimes  can  be  expressed  in  the  form  of 
a  drag  force: 

Fr  -  (l/2>oeC0Af  (Vf  -  v/  (4-20) 

where  pg  =  mass  density  of  the  earth  cover 

Cq  =  drag  coefficient 

=  effective  cross-sectional  area  of  the  fragment 

1/^  =  speed  of  fragment 

V  =  mean  speed  of  cover 
w 

There  are  at  least  two  convenient  ways  to  approximate  the  motion  of  the 

earth  cover  after  detonation.  One  method  is  to  assume  that  the  thickness 

of  the  earth  cover  remains  constant  as  the  cover  moves;  this  specifies 

2 

that  the  density  of  the  cover  decreases  as  (1/R  )  in  order  to  conserve  mass. 
The  other  view  is  to  assume  that  the  density  of  the  cover  is  a  constant, 
so  that  the  thickness  of  the  cover  decreases  as  (1/R  )  as  the  cover  moves  out. 

Solutions  of  the  equations  of  the  pre-interaction  regime  specifies  at 
what  radius  and  time  the  interaction  between  a  fragment  of  mass  m  and  the 
earth  cover  will  commence;  the  solution  yields  the  initial  values  required 
for  the  solution  of  the  equations  of  the  interaction  regime.  The  inter¬ 
action  equation  determines  the  velocity  with  which  the  fragment  leaves 
the  wedge  after  passing  through  it.  It  should  be  obvious  that  some  fragments 
will  not  penetrate  through  the  cover.  Ry  considering  the  masses  and  velo¬ 
cities  of  fragments  as  specified  by  these  distributions  for  the  weapon,  it 
is  possible  to  estimate  the  fragments  which  will  and  will  not  penetrate  the 
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cover.  It  can  be  seen  that  the  effect  of  the  cover  is  to  alter  the  mass  and 
velocity  distributions  of  the  fragments.  These  new  distributions  are  assumed 
again  to  be  the  initial  ones  for  a  single  equivalent  weapon,  and  estimations 
of  far-field  fragment  densities  and  damage  levels  proceed  according  to  GARD's 
current  technique. 

Based  on  these  comments,  the  first  model  to  approximate  the  effect 
of  the  earth  cover  is  presented.  The  earth  cover  moves  out  very  slowly 
compared  with  the  fragments  from  the  stack  and  the  fragments  from  the  stack 
must  perforate  an  earth  layer  which  is  essentially  intact.  The  earth  cover, 
is  however  broken  up  and  behaves  essentially  as  a  fluid.  The  main  effect 
of  the  cover  is  increased  drag  on  the  fragments.  Experiments  by  Allen  et 
al.  (1957a,  1957b)  on  the  retardation  of  steel  spheres  by  sand  show  that, 
for  high  fragment  velocities,  the  main  retardation  is  essentially  of 
this  form;  i.e.,  the  resisting  force  on  a  fragment  is 

F  =  (7)  p  CD  V2A  =  B*m1/3V2 

where  p$  is  the  sand  density 

CQ  is  1.6  for  spheres  perforating  sand 

V  is  the  fragment  velocity 

A  is  the  presented  are  of  the  fragment 

6*  =  Bp  /p  • 

‘V^air 

We  assume  that  the  soil  cover  resists  the  fragment  in  the  same  way 
as  sand.  In  this  model,  the  effect  of  the  cover  is  to  reduce  the  fragment's 
velocity  as  it  passes  through  the  cover.  Since  the  cover  does  not  move 
appreciably  during  the  perforation  process,  we  assume  it  is  stationary  and 
that  the  fragment  must  penetrate  the  original  thickness.  Further  we 
can  neglect  the  cover  velocity  in  calculating  the  resisting  force  on  a 


GFNERAL  AMERICAN  RESEARCH  DIVISION 

60 


fragment.  Denoting  the  distance  from  the  munition  stack  to  the  inner  wall 
of  the  igloo  by  d,  and  the  thickness  of  the  igloo  by  h  (where  d  and  h  are 
functions  of  the  initial  angle  and  azimuth  angle)  the  fragment  velocity  at 
the  time  the  fragment  strikes  the  inner  wall,  is 


V1  =  VQ  exp  (-gm“1//3d) 


■0  ^  =  VQ  exp  (-e*m'1/3h) 


(4-22) 

The  fragment  velocity  when  it  emerges  from  the  earth  layer  is 

V2  =  V]  exp  (-e*m“1/3h)  (4-23) 

By  comparing  the  fragment  velocity,  V^,  with  the  fragment  velocity 
at  a  distance,  r  +  h,  in  air  alone,  we  define  an  effective  initial  velocity 

V 

The  effect  of  the  earth  cover,  in  this  model  is  to  retard  the  fragment; 

the  magnitude  of  this  effect  can  be  represented  as  a  fictitious  lower 

initial  velocity  on  the  unit  hemisphere.  Using  the  value  of  Pearth  “ 

3 

100  lb/ft  ,  Cg  =  1.60  and  the  value  of  h  as  a  function  of  elevation  angle 
from  the  center  of  the  igloo,  we  can  readily  calculate  the  far-field  fragment 
density  (we  assume  the  hemisphere  data  for  the  side-spray  direction). 

Figure  4-8  through  4-11  show  the  reduction  in  far-field  fragment  density 
for  a  single  munition.  This  model  is  denoted  as  the  penetration  model  on 
the  figures. 

There  is  yet  another  way  to  consider  the  action  of  the  earth  cover  as 
it  moves  out;  i.e.,  it  breaks  up  into  chunks  rather  than  infinitesimal 
particles.  The  problem  then  becomes  very  similar  to  that  of  the  two  side- 
by-side  munitions  whose  fragments  interact.  This  approach  involves  calcu¬ 
lation  of  the  probability  of  collision,  and  the  effects  of  the  interaction 
process  are  calculated  essentially  identically  to  the  method  described  in 
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Fiqurt  4-8,  EFFECT  OF  EARTH  COVER  ON  FRAGMENT  DENSITY  VERSUS  RANGE 
MI07  155mm  PROJECTILE  TNT  LOADED  (SIDE  DIRECTION) 


FRAGMENT  DENSITY 


Figure  4-  9  ,  EFFECT  OF  EARTH  COVER  ON  FRAGMENT  DENSITY  VERSUS  RANGE 
M437A2  175mm  PROJECTILE  (SIDE  DIRECTION) 


FRAGMENT  DENSITY 


RANGE  -  k  ft. 


Figur#  4-10,  EFFECT  OF  EARTH  COVER  ON  FRAGMENT  DENSITY  VERSUS  RANGE 
MK  82  500-1*  r  'JMB  (SIDE  DIRECTION) 


GM  ENT  DENSIT  Y 


the  preceding  section.  Results  using  this  method  are  senitive  to  the 
mass  distribution  assumed  for  the  chunks  of  the  earth  cover.  The  mathematical 
models  developed  for  fragment-fragment  interaction  can  be  used  for  the 
fragment-cover  interaction. 

This  second  model  again  assumes  the  cover  doesn't  move  very  far  when 
the  fragments  overtake  it,  but  the  interaction  of  fragments  with  the  cover 
is  described  differently.  The  cover  is  assumed  to  break  up  into  fragments 
that  have  a  mass  distribution  similar  to  the  fragment-distribution.  If 
a  bomb  fragment  strikes  an  earth  clod  which  is  larger,  it  is  stopped 
completely.  If  it  strikes  a  cold  which  is  smaller,  it  is  not  affected. 

This  model  is  a  simplification  of  completely  inelastic  impact  in  which  the 
velocity  of  the  fragment  and  clod,  upon  impact,  have  the  same  final 
velocity. 


VfMf  +  VeMe 


where  Vf  is  the  final  velocity  of  both 

Vp  is  the  fragment  velocity  at  impact 

Vg  is  the  earth  clod  velocity  at  impact 

Mp  is  the  mass  of  the  fragment 

is  the  mass  of  the  earth  clod 
e 

The  effect  of  this  model  is  to  reduce  the  number  of  fragments  (on  the 

unit  hemisphere)that  are  thrown  out.  Mathematically  this  is  represented 

by  an  altered  mass  distribution  f.(m) 

1  1/3 

f^m)  =  f(m)[l  -  exp  (-(™)  yh)] 

where  y  represents  the  probability  of  impact  of  a  fragment  by  an  earth  clod 


when  the  fragment  traverses  one  foot  of  earth. 
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It  is  very  hard  to  estimate  accurately  the  value  of  y  to  be  used  in 
this  equation.  We  take  it  as  the  probability  that  a  metal  fragment  will 
strike  at  least  one  earth  fragment  in  j  foot  penetration  of  the  earth  cover. 

Also  we  assume,  for  lack  of  better  estimates,  that  the  u  of  the  earth  clod 

distribution  is  the  same  as  that  for  the  fragment  distribution.  With  these 

crude  estimates  we  show  the  far-field  fragment  densities  for  this  type  model 

in  Figures  4-8  through  4-11  for  a  single  munition  and  the  standard  igloo. 

This  model,  with  the  very  optimistic  values  for  a  fragment  getting 
through  without  impact,  gives  very  low  values  for  the  far-field  fragment 
density. 

Under  the  effect  of  the  blast,  the  cover  breaks  up  and  moves  out. 

Thus  a  third  type  of  interaction  must  be  considered  wherein  the  cover  frag¬ 
ments  move  essentially  at  the  same  velocity  as  the  munition  fragments  and 
the  number  of  collisions  is  minimal.  The  third  model  differs  from  the  other 
two  in  that  it  is  assumed  that  the  earth  cover  breaks  up  and  moves  away  fast 
enough  so  that  no  interaction  between  metal  fragments  and  the  cover  will 
occur  for  any  initial  elevation  angle  greater  than  the  final  crater  lip. 

For  elevation  angles  lower  than  the  final  crater,  the  blast  is  trying  to 
move  a  very  large  amount  of  earth  and  will  not  break  it  up.  For  the  higher 
angles,  the  cover  is  assumed  completely  ineffective  in  stopping  fragments 
and  for  the  lower  angles,  all  fragments  are  considered  to  be  stopped. 

Figure  4-12  schematically  illustrates  this  model.  Estimates  of  a  ,  the 

m 

angle  to  the  top  of  the  residual  cover  are  obtained  from  the  profiles  of 

the  craters  in  the  Arco  detonations  and  Eskimo  I.  An  estimate  of  a  'v  10° 

m 

is  consistent  with  this  data.  In  this  model  then,  the  calculation  of  the 
fragment  density  in  the  far-field  is  effected  by  replacing  the  lower  limit 
in  the  integration  over  a0  in  Equations  3-21  or  3-22  by  am» 
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Figure  4*12,  SCHEMATIC  DIAGRAM  ILLUSTRATING  THE  THIRD  MODEL 
OF  EARTH  COVER  EFFECT  ON  FRAGMENTS- 
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4.4  Other  Parameters 


The  effect  on  the  far-field  fragment  density  when  two  of  the  initial 
fragment  parameters  are  altered  is  now  examined.  As  noted,  the  fragment 
mass  distribution  seems  to  be  changed  «nen  a  stack  of  munitions  is  detonated 
with  a  single  munition  primed  (Draper  and  Watson,  1970,  Feinstein  and 
Nagaoka,  1970,  Feinstein,  1972b).  The  observations  indicate  that  more  large 
fragments  are  generated.  There  is  no  available  data,  however,  on  the  mass 
distribution  that  is  obtained;  therefore,  a  simple  model  to  examine  this 
effect  is  used.  Assume  the  altered  mass  distribution  has  the  same 
functional  form  as  the  arena  data  but  with  fewer  total  fragments  and  a 
larger  mean  fragment  mass.  This  is  effected  in  the  formula  by  changing 
the  value  of  u  to  yu,  where  y  is  a  constant  factor  and  by  changing  N  to 

No  0 

—  ,  to  conserve  the  total  mass  of  fragments  emitted.  The  effect  of  this 
change  is  illustrated  for  fragment  densities  in  the  side  direction  for  the 
Ml  1 7  750  lb  bomb.  Figure  4-13  shows  the  fragment  density  versus  y  for 
1  •<  y  <_  10  at  selected  ranges  between  1000  ft  and  2500  ft.  For  the  smaller 
range,  the  fragment  density  decreases  with  increasing  y,  decreasing  about 
an  order  of  magnitude  as  y  increase  about  an  order  of  magnitude.  At  the 
other  extreme  (i.e.,  ranges  of  2000  and  2500  feet),  the  change  in  fragment 
density  with  y  is  very  slight,  being  about  30%  over  the  range  of  y  presented. 
Although  there  are  fewer  total  fragments  for  larger  y,  there  is  a  larger 
fraction  of  large  fragments  and  the  larger  fragments  are  hurled  to  larger 
distances  than  smaller  fragments. 

As  far  as  the  far-field  fragment  hazard  is  concerned,  the  effect  of 
changing  the  mass  distribution  in  a  stack  detonation  is  slight. 

Now  turn  to  a  feature  whose  effect  will  not  be  slight  in  estimating 
the  far-field  fragment  hazard.  In  solving  the  ballistic  trajectory  equations, 
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it  is  normally  assumed  (in  all  the  theoretical  works  presented)  that  the 
fragments  are  spinning  as  they  leave  the  detonation  site  and  that  the 
presented  area  can  be  represented  (at  least  over  the  larger  trajectories) 
as  the  mean  presented  area.  It  would  be  possible  for  some  fragments  to  be 
ejected  with  no  or  little  spin  or,  if  spinning  at  the  outset,  to  stabilize 
with  their  major  axis  roughly  parallel  to  the  trajectory.  In  this  case  the 
presented  area  would  be  minimized  and  the  ballistic  density  maximized. 

For  example,  look  at  a  fragment  with  square  cross-section  and  with  an 

aspect  ratio  of  2:1,  aligned  with  the  trajectory.  For  steel  fragments,  the 

3 

ballistic  density  is  3960  gr/in  or  six  times  the  mean  value. 

The  far-field  fragment  density  was  calculated  for  various  values  of 
the  ballistic  density,  ranging  from  the  mean  value  given  in  the  arena  data 
to  four  times  the  arena  data  (i.e.,  from  ^  600  gr/in  to  2500  gr/in  ). 

Figure  4-14  shows  the  fragment  density  versus  ballistic  density  at  selected 
ranges  from  1000  feet  to  2500  feet  for  the  Ml  17  750  lb  bomb  in  the  side 
spray  direction.  At  the  larger  ranges  (e.g.  2500  feet)  a  dramatic  increase 
in  the  fragment  density  is  noted,  increasing  by  almost  two  orders  of  magni¬ 
tude  as  the  ballistic  density  changes  by  a  factor  of  four.  It  is  suggested 
that  this  factor  is  quite  important  for  far-field  fragment  estimations. 

Now  it  is  improbable  that  all  fragments  will  be  stabilized  and  oriented  to 
the  trajectory,  but  even  if  only  10%  of  them  are,  while  the  remaining 
90%  are  spinning;  the  increase  of  the  far-field  fragment  density  will  be 
dramatic.  This  increase  in  effective  ballistic  density  could  also  be 
interpreted  as  a  smaller  effective  drag  coefficient  on  the  fragment  during 
its  travel.  Increases  in  the  drag  coefficient  by  50%,  equivalent  to  changes 
in  the  ballistic  density  by  a  factor  of  3,  are  unlikely. 
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The  effect  of  increased  initial  velocity  observed  by  Draper  and 
Watson  (1970)  is  less  dramatic  in  increasing  the  far-field  fragment  density. 
An  examination  of  the  maximum  range  of  a  fragment  of  mass  explains  this. 

The  maximum  range  is 

.2/3,.,  1/3  1/6 

R  =  ti(~)  m  /J[.5140  -  1.0368  ln(137(— )  )]  (4.24) 

KA  kA  Vo 

where  k  is  the  actual  ballistic  density 

is  the  arena  data  average  ballistic  density 
VQ  is  the  initial  velocity 

The  dominant  dependence  on  VQ  is  with  In  VQ  while  the  dominant  dependence 
2/3 

on  <  is  as  <  . 

Note  changes  in  VQ  by  a  factor  of  2  or  3  may  be  expected,  while  the 
value  of  k  for  selected  projectiles  could  range  easily  up  to  5-10  times 
the  arena  value. 
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Chapter  V 


COMPARISON  OF  THE  THEORETICAL  MODELS  WITH  EXPERIMENTAL  DATA 


5.1  Stack  Models  and  Measurements 


The  elementary  theory  suggests  that  the  effective  number  of  munitions 
is  related  to  the  number  of  munitions  in  the  surface  layers  on  the  top  and 
the  exposed  side  where 

N  ff  =  0.9NS  +0.1  Nt 

First,  the  theoretical  calculations  are  compared  with  the  fragment  density 
in  the  far-field  (500  feet  to  1500  feet)  for  the  side  spray  direction  from 
tne  2x5  and  5x3  stacks  of  Ml  17  750  lb  bombs  tested  at  NWC  (Feinstein  and 
Nagaoka,  1970).  In  the  former  N-j.  =  2,  and  Ny  =  5  in  the  latter,  while 
N<.  =  3  for  both.  Thus 


j  2.9  for  2  x  3  stack 
Neff  *  \ 

eTT  \  3.2  for  5  x  3  stack 

Figure  5-1  shows  the  comparison  using  the  arena  data  approximation.  The 
comparison  is  quite  good. 

The  next  comparison  is  with  the  first  two  events  in  the  Big  Papa  test. 

These  tests  were  the  detonation  in  an  open  revetment  of  a  mixed  stack  of 

Ml  1 7  750  lb  bombs  and  M66A2  2000  lb  bombs.  The  750  lb  and  2000  lb  bomb  are 

both  tritonal  filled  and  have  the  same  C/M  ratio.  The  arena  data  for  the 

2000  lb  bomb  has  the  same  fragment  velocity  (within  10%)  and  the  same 

8 

average  fragment  mass  as  the  750  lb  bomb  with  approximately  y  more  fragments. 
Therefore,  each  2000  lb  bomb  is  replaced  by  2-2/3  750  lbs  for  the  purpose 
of  comparison. 
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RANGE 


Figur*  5-1,  COMPARISON  OF  CALCULATED  FRAGMENT  DENSITY  VERSUS 

RANGE  WITH  EX  PE  RIM  ENT- SMALL  STACK  OF  Mll7  750-lb  BOMBS 
(SIDE  SPRAY  DIRECTION) 


The  stacks  for  these  detonations  were  rectangul ar.  Details  of  the 
stack  composition  are  given  in  Peterson  (1968).  In  the  first  test,  there 
were  64  750-lb  bombs  in  the  top  l^yer  and  4  750-lb  and  12  2000-lb  bombs 
in  the  side.  For  the  second  test,  there  were  21  750-lb  bombs  in  the  top 
and  4  750-lb  and  8  2000 -T b  bombs  on  tne  side. 

Thus  =  39  in  the  first  case  and  25  in  the  second.  The  comparison 

between  theory  (using  arena  data)  and  measured  values  is  not  very  good; 
between  2000  ft  and  3000  ft  there  is  roughly  a  factor  of  50  between  the 
two.  Schreyer  and  Romesberg  (1970)  noted  that  using  arena  data  and  a 
multiplier,  an  effective  stack  of  267  bombs  was  required  to  match  the  data 
between  2000  and  3000  ft. 

Peterson  (1968)  measured  the  mass  of  the  smallest  fragments  at  1800  and 

3000  ft.  These  fragments  are  much  smaller  than  that  predicted  from  the 

ballistic  equations.  The  initial  velocity  required  for  these  fragments  is 

well  in  excess  of  25,000  ft/sec.  However,  if  the  effective  ballistic  density 

was  twi'e  the  arena  average  these  small  fragments  could  reach  these  distances 

with  no  increase  in  initial  velocity.  Therefore,  the  fragment  densities, 

•  3 

were  recalculated  using  N  ^  and  k  =  1200  grains/in  .  Figure  5-2  shows  the 
comparison  of  the  predicted  and  experimental.  The  comparison  is  not  great 
(a  factor  of  'v  2  at  2000  ft),  but  probably  reasonable.  There  are  very 
slight  changes  in  q  from  500  to  1500  ft.  suggesting  that  the  agreement 
with  the  NWC  2x3,  5x3  tests  remains  valid. 

A  more  realistic  source  model  would  probably  have  the  smaller  fragments 
starting  out  with  a  higher  velocity  than  the  larger  fragments  and  would  also 
exhibit  dependence  of  the  shape  factor  on  the  size.  Since  there  is  no 
data  available  on  this  dependence,  no  attempt  to  fit  any  such  variation  was 
made. 
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FRAGMENT  DENSITY 


5.2  Igloo  Models  and  Measurements 

There  is  only  one  experiment  for  weapons  contained  within  an  igloo  that 
can  be  compared  with  the  theoretical  far-field  fragment  densities.  The 
Eskimo  I  test  of  13696  155mm  projectiles  detonated  in  a  standard  Army  Igloo 
were  measured.  The  projectiles  were  placed  in  a  stack  shown  in  Figure  5-3. 
The  munitions  were  oriented  vertically.  Figure  2-5  shows  the  plan  view  of 
the  test,  and  indicates  the  lines  along  which  fragment  measurments  were  made. 
The  number  of  the  projectiles  on  the  top  layer  for  considering  fragments 

off  the  end  is  3700;  off  the  side  almost  half  of  these  are  screened  by  the 
odd  shape  of  the  stack;  thus  =  2300  there.  Off  the  end  Ns  =  240  and 

off  the  side  N<.  =  230.  First  the  far-field  fragment  densities  from  this 
stack  were  calculated  as  if  the  stack  were  in  the  open,  using  the  isotropic¬ 
ally  averaged  arena  data.  These  values  are  quite  low  beyond  1000  ft.  This 
predicted  density  has  a  Sq/aR  greater  than  that  measured.  The  two  models 
with  strong  cover- fragment  interaction  will  have  an  even  greater  3q/3R 
and  would  thus  deviate  further. 

The  mass  distribution  of  the  fragments  were  measured  at  each  range. 
(Weals,  1973,  Feinstein,  1972b).  Figure  5-4  shows  the  minimum  mass  fragment 
recovered  as  a  function  of  range  for  this  experiment.  Also  plotted  is 

3 

the  theoretical  minimum  mass  versus  range  using  k  =  660  grams/in  , 

VQ  =  5000  fps;  k  =  1200  grains/in3,  VQ  =  5000  fps;  and  k  =  660  grains/in3, 

VQ  =  40,000  ft/sec.  In  line  with  the  previous  fit  to  the  Big  Papa  experiment 
the  velocity  increase  to  account  for  the  mass  is  rejected  as  the  prime 
explanation  and  the  ballistic  density  cor ''espon ding  to  the  more  slender 
projectile  is  accepted. 

3 

The  calculations  were  repeated,  but  a  value  of  k  *  1200  gr/in 
was  used.  The  calculated  values  agree  to  within  a  factor  of  two  for  most 
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Figurt  5- 


CROSS  SECTION  OF  PALLET  STACK  IN  IGLOO  MAGAZINE 
ESKIMO  I. 


9 


Figure  5-4,  MEASURED  MINIMUM  MASS  FRAGMENTS  VERSUS  RANGE 
ESKIMO  I  (EXP.  WEALS  1973)  COMPARED  WITH  THEORY-’ 
FRAGMENTS  EJECTED  AT  OPTIMUM  LAUNCH  ANGLE. 
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of  the  range  and  the  slope  is  qualitatively  in  agreement.  The  perforation 
model  for  the  igloo  effect  would  yield  values  lower  than  these  by  better 
than  an  order  of  magnitude.  The  model  wherein  fragments  are  stopped  for 
low  elevation  angles,  but  otherwise  unimpeded,  yields  values  which  are 
in  good  agreement  with  the  measurements.  Figure  5-5  shows  the  measured 
data  and  the  theoretical  predictions  for  this  model.  (For  fragments  in 
the  direction  of  the  headwall ,  assume  the  headwall  is  destroyed  by  the 

blast  and  offers  no  resistance  to  low  elevation  angle  fragments  in  that 

direction.)  For  this  model  the  N  ^  are  as  follows: 

Headwall  N  ff  =  0.9  (240)  +  0.1  (3700)  =  586  *  600 

Side  N  ff  =  0.7  (280)  +  0.1  (2300)  =  426  *  425 

Back  N  ff  =  0.7  (240)  +  0.1  (3700)  =  538  550 

5.3  Comments 

The  theoretical  predictions  for  far-field  fragment  densities  for  the 
accidental  detonation  of  stacked  munitions  and  of  stacked  munitions  within 
igloo  magazines  compares  quite  well  with  the  few  adequate  measurements  of 
the  far-field  fragment  densities. 

A  commentary  of  the  assumptions  that  are  used  and  of  the  models  that 
were  developed  that  best  fit  the  data  is  useful. 

i.  The  fragment  densities  in  the  far-field  are  calculated  using 
an  approximation  :o  the  exact  solution  to  the  ballistic  trajectory  equations. 
The  approximation  technique  is  very  good  for  the  fragment  masses  and  initial 
velocities  of  interest.  The  ballistic  equations  that  this  approximation 
represents  ignores  any  yaw  effects  and  any  wind  conditions.  The  spin  of 
the  fragments  is  approximated  by  assuming  an  effective  mean  presented  area 
only  and  no  deviations  from  the  spinless  trajectory  by  Magnus  and  cross  forces 
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is  considered.  This  effect  of  spin  on  the  trajectory  was  not  estimated, 
but  it  could  be  appreciable. 

2.  The  ballistic  equations  allow  the  fragment  density  to  be  computed 

at  large  ranges  from  the  specification  of  fragment  number,  mass,  initial  velo¬ 
city  and  ballistic  density  at  any  specified  azimuth  and  elevation  on  the  unit 
hemisphere  centered  at  ground  zero.  Given  a  set  of  arena  data  for  a  single 
munition,  the  far- field  fragment  densities  may  be  readily  computed. 

3.  The  model  for  the  stacks  led  to  a  definition  of  an  effective 
number  of  single  munitions.  Only  munitions  in  the  top  layer  and  the  side 
layer  nearest  the  observation  point  contribute  significantly  to  the  far- 
field  fragment  density.  The  bulk  of  the  contribution  comes  from  the  side 
layer.  Fragments  from  munitions  in  the  interior  of  the  stack  are  effec¬ 
tively  stopped  by  the  fragments,  blast,  and  combustion  products  from  the 
adjacent  munitions. 

4.  The  best  model  for  describing  the  effect  of  an  earth  cover  on 

the  retardation  of  fragments  is  that  the  crown  of  the  earth  cover  is  blown 
apart  by  the  blast  and  moves  out  as  part  of  the  fragment  pattern  and  does 
not  significantly  interact  with  the  munition  fragments.  At  very  low 
elevation  angles  the  cover  essentially  remains  in  place  and  completely 
stops  all  fragments  with  very  low  elevation  angles.  The  original  purpose 
for  the  design  of  the  earth  cover  is  the  prevention  of  intermagazine 
communication  of  the  detonation  by  fragment  impact  and  blast.  As  far  as 
fragments  are  concerned,  this  objective  is  met.  The  cover  affords  little 
reduction  in  the  far- field  fragment  numbers,  though. 

5.  The  distribution  of  fragment  masses,  total  fragment  numbers  and 
fragment  initial  velocities  from  an  accidental  detonation  in  a  stack,  where¬ 
in  one  of  the  munitions  detonates  and  the  remainder  are  detonated 
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sympathetically,  may  differ  significantly  from  the  values  obtained  by 
detonating  a  single  munition.  In  particular,  larger  fragments  and  higher 
initial  velocities  may  be  expected.  The  effect  of  larger  fragments  (to¬ 
gether  with  its  corollary  property  of  fewer  of  them)  do  not  significantly 
alter  the  densities  that  may  be  expected  in  the  ranges  1000-3000  feet. 
Increases  in  initial  velocity  affect  the  fragment  densities  only  slightly. 
The  parametric  studies  indicate  that  the  fragment  densities  in  the  far-field 
are  extremely  sensitive  to  effective  ballistic  density.  A  change  in  this 
variable  by  a  factor  of  two  leads  to  increases  in  the  number  of  fragments 
observed  in  the  far-field  by  better  than  an  order  of  magnitude. 

6.  Use  of  arei.a  data  could  not  effect  a  match  between  theory  and 
the  measured  data  from  the  Big  Papa  and  Eskimo  I  experiments.  The  de¬ 
tailed  measurements  of  these  experiments  show  that  much  smaller  fragments 
are  ejected  to  larger  distances  than  could  be  expected  using  the  initial 
velocities  and  the  average  shape  factor  from  the  arena.  Modifications  of 
either  of  these  factors  could  account  for  the  observed  data.  The  velocity 
increments  required  are  large,  being  a  factor  of  3-10  times  the  maximum 
arena  velocity,  while  only  small  changes  in  the  average  effective  shape 
factor  are  required. 

7.  An  additional  assumption  for  use  in  the  theoretical  calculations 
at  this  stage  of  development  is  that  the  initial  hemisphere  parameters 
should  be  chosen  to  be  isotropic  with  the  largest  values  of  fragment 
number  and  velocity  from  the  arena  data  chosen.  The  munitions  in  an 
accidental  detonation  will  be  tossed  about  prior  to  initiation  of  the 
sympathetic  detonation  and  the  actual  orientation  will  be  subject  to  some 
sort  of  statistical  distribution.  This  project  is  interested  in  the 
determination  of  safe  distances  and  to  be  conservative  it  should  be 
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assumed  that  the  maximum  source  is  operative  in  any  direction. 

It  is  worth  noting  that  either  greatly  enhanced  initial  velocity  or 
an  altered  ballistic  density  is  required  to  explain  the  measured  far-field 
fragment  density.  We  lean  toward  the  latter  as  the  fragment  density  is 
extremely  sensitive  to  small  changes  in  the  shape  factor,  while  extremely 
large  initial  velocity  increments  are  required.  If  the  latter  case  is 
descriptive,  some  thought  must  be  given  to  the  injury  criteria.  Fugelso 
(Fugelso,  et  al.  1960,  1961,  1966)  has  shown  that  perforation  limit  velo¬ 
cities  for  thin  metallic  plates  are  highly  dependent  on  the  mass  to  pre¬ 
sented  area  ratio  at  impact.  Sperrazza  and  Kokinakis  (1968),  and  Kokinakis, 
(1970)  have  shown  similar  dependence  for  perforation  limit  velocities  of 
skin.  Their  data  yield 

a 

V1 imi t  “  a  +  ^  k  a»  b»  constants 

The  trauma  data  (Bowen  et.  al .  1968)  for  unilateral  and  bilateral 

lung  hemorrhage  by  blunt  missile  impact  is  based  on  a  set  of  data  with 
two  sizes  of  projectiles  ,M  ■  .4  and  .8  lb, with  A  constant  for  both.  The 
skull  fracture  data  (Feinstein,  et  al.  1968,  and  Feinstein  1971),  is 
based  on  dropping  skulls  on  a  flat  surface,  not  by  impacting  them  with  a 
projectile  at  all  (Gurdjian  1949). 

In  view  of  the  high  possibility  of  the  far-field  fragments  being  highly 
elongated,  it  is  suggested  that  the  critical  velocity  requirements  be 
reexamined  quite  thoroughly.  For  tr.e  present  it  is  suggested  that  all 
fragments  beyond  1000  feet  be  treated  as  potentially  lethal. 
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Chapter  VI 

CONCLUSIONS  AND  CONJECTURES 


6.1  Summary  of  the  Model  for  Far-Field  Fragment  Density  Calculations 

On  the  basis  of  the  theory  and  their  limited  comparisons  with  experi¬ 
mental  data,  the  following  model  for  the  calculation  of  the  far-field  frag¬ 
ment  density  from  the  accidental  detonation  of  stacks  of  munitions  is 
tentatively  proposed. 

I.  Each  munition  generates  an  isotropic  fragment  distribution  with 
the  total  number  of  fragments  from  the  arena  data  distributed  uniformly 
over  a  unit  sphere.  These  numbers  have  an  upper  bound  given  by  the  average 
densities  in  the  arena  data  in  the  side  direction.  Since  the  calculations 
are  to  be  conservative  for  safety  purposes,  we  use  these  av'  .  ;es  (Table  2). 
The  average  mass  will  be  taken  as  that  of  the  arena  data.  The  initial 
fragment  velocity  will  be  the  maximum  verity  of  the  arena  data,  this  is 
normally  the  maximum  side-spray  velocity.  The  ballistic  density  will  be 
taken  as  twice  that  given  in  the  arena  data. 

II.  Let  Ny  denote  the  number  of  munitions  in  the  top  layer  of  the 
stack  and  N„  the  number  of  munitions  in  nearest  side.  The  total  number 

o 

of  effective  munitions  is: 

1)  For  an  open  stack 

Neff  =  0,9  NS  +  0,1  NT  (6-1) 

2)  For  a  stack  within  a  standard  earth-covered  igloo 

Neff  =  0,7  NS  +  0,1  NT  (6-2) 

In  the  headwall  direction,  treat  it  as  in  the  open. 
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III.  With  these  numbers  as  input,  the  far-field  fragment  density  is 
calculated  by  the  approximate  formula 

^stack  ~  ^eff  ^single 
and  qsingie  is  given  by  Equation  (3-22). 

IV.  All  fragments  are  considered  hazardous.  Figure  6-1  through 
6-4  presents  the  fragment  densities  for  the  single  munition  for  the 
four  munition  types  considered. 

6.2  Tentative  Quantity-Distance  Relationships 

Tentative  quantity-distance  curves  for  the  four  munitions  are 
presented.  An  assumption  on  typical  stack  shapes  are  required  to  illus¬ 
trate  the  calculations.  Let  M  equal  the  total  number  of  munitions  in  the 
rectangular  stack  with  dimensions  3n  x  n  x  n  (Figure  6-5).  For  this 
example 

Nt  =  3^ 

and  Nj.  =  3n2  off  the 

=  n2  off  the 

Thus  =  3n2  off  the 

=  1.2  n2  off  the 

=  2.4  n2  off  the 

=  n2  off  the 

and  M  =  3n3 

The  quantity  of  explosive,  W,  is 
W  =  MW] 

where  is  the  explosive  weight  in  a 


end 

side 

side 

>  open  store 

end 

/ 

si  de 

\  stored  in  earth- 

end 

/  covered  igloo 

single  munition. 
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FRAGMENT  DENSITY 


Figures  6-6  through  6-9  show  the  quantity-distance  relations  for  these 
four  cases  for  the  four  munitions.  Values  of  the  fragment  density  of 
1/600  square  feet  and  1/6000  square  feet  are  shown  on  each  graph.  Also 
plotted  on  each  curve  a  a  the  inhabited  building  distarce  for  Class  7 
explosives  (D0D  Manual  4145,  27M,  1971)  and  the  British  criteria  R  =  515  W^5 
for  a  10~5  probability  of  being  struck  by  a  fragment  (Jarrett  1968).  It 
should  be  noted  that  the  tentative  quantity-distances  curves  for  the 
fragment  hazard  indicate  that  the  fragment  hazard  might  be  the  controlling 
safety  feature  for  stacks  with  less  than  100,000  to  200,000  lbs  of  explosive. 
The  quantity- distance  relation  for  these  fragment  hazards  will  increase  the 
required  safe  distances  for  the  smaller  stacks. 
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EXPLOSIVE  WEIGHT, - - tb*. 


Figure  6*6,  QUANTITY*  DISTANCE  FOR  STACKS  OF  MI07  155  mm  PROJECTILES 


EXPLOSIVE  WEIGHT 


Figure  6-8,  QUANTITY-DISTANCE  FOR  STACKS  OF  MK82  5C0-lb. BOMBS. 
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* 


Figure  6-9,  QUANTITY- DISTANCE  FOR  STACKS  OF  M  1 17  780-lb.  BOMBS 


Cnapter  VII 
RECOMMENDATIONS 


During  the  course  of  this  v.'ork,  the  lack  of  a  consistent  set  of 
experimental  data  on  fragments  from  large  size  stacks,  both  contained  and 
in  the  open,  hampered  conclusive  evaluation  of  the  theoretical  work.  This 
lack  prohibits  any  generalizations  or  extrapolations  to  be  made  with  any 
degree  of  confidence.  The  need  for  a  well  planned,  consistent  experimental 
program  for  fragment  hazards  from  munition  stacks  is  paramount  if  any 
form  of  quantitative  predictions  are  ever  to  be  made. 

The  original  intent  of  this  work  was  to  define  the  quantity-distance 
relationship  for  fragment  hazards  in  terms  of  the  munitions,  stack  and 
cover  parameters  in  such  a  manner  that  this  information  could  be  displayed 
on  slide  rule  or  nomograph  in  a  manner  similar  to  the  quantity-distance 
relationships  for  blast  effects  from  accidental  detonation  of  munition 
stores  (Fugelso  et  al  •  1972a, b).  This  current  work  has  progressed  to  a 
point  where  the  parameters  have  been  identified  and  good  theoretical 
agreement  has  been  obtained  with  a  very  limited  number  of  tests.  While 
the  agreement  to  date  is  encouraging,  it  is  not  conclusive  enough  to  commit 
it  to  a  general  computation  aid. 

To  start,  a  consistent  experimental  program  for  the  fragment,  patterns 
from  munitions  in  open  stores  should  be  made.  The  program  consists  of  many 
detonations  of  a  single  type  of  munition,  in  a  variety  of  configurations. 
First,  the  far-field  fragment  distribution  from  a  single  munition  should 
be  made,  in  both  the  horizontal  and  vertical  orientations.  This  plan  is 
necessary  to  determine  the  effect  of  initial  fragment  shape 
distribution.  (Remember,  the  trajectory  equations  act  as  high  ballistic 
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density  and  high  initial  velocity  filters  for  long  range  propagation  and  the 
distances  traversed  by  a  fragment  are  extremely  sensitive  to  small  variations 
in  these  parameters.  The  arena  data  as  presented,  gives  only  average 
values  for  these  parameters.)  Sufficient  duplication  of  the  experiment  must 
be  made  to  give  statistical  significance  to  the  far-field  data. 

From  the  tests  of  the  individual  munitions,  we  progress  to  stacks.  A 
sufficient  number  of  stack  tests  must  be  made  to  determine  the  surface, 
bulk  and  orientation  dependence  on  long  range  propagation.  Increased 
initial  velocities,  altered  fragment  distributions  and  the  tendency  to 
source  isotropy  must  be  evaluated.  It  is  obvious  that,  in  addition  to 
far-field  fragment  measurements,  measurements  of  the  distribution  of 
initial  velocity  and  fragment  distribution  at  the  source  must  be  made. 

As  a  part  of  the  stack  tests,  tne  effects  of  inter-munition  spacing, 
munition  orientation,  and  mode  of  initiation  of  the  stack  (e.g.,  all  munitions 
primal,  or  the  munition  primed  In  the  center  or  at  the  edge  of  the  stack) 
must  also  be  evaluated.  Again,  sufficient  duplication  of  these  experiments 
Is  essential. 

Upon  conclusion  of  these  tests,  the  essential  parameters  of  the  sources 
as  a  function  of  stack  parameters  must  b;?  determined  and  the  dependence 
of  the  far-field  fragnent  densities  as  a  function  of  these  must  be  demon¬ 
strated. 

Then,  and  only  then,  can  the  effect  of  the  earth  cover  be  tested,  In 
the  same  manner.  First,  munition  in  revetments  should  be  tested,  followed 
by  larger  stacks  In  Igloos.  It  Is  necessary  In  the  Igloo  test  that  the 
destruction  of  the  Igloo  be  accomplished  in  the  same  manner  as  it  would  be 
In  a  larger  scale  accidental  detonation. 
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Another  area  for  possible  study  is  the  criterion  for  damaging  fragments 
when  the  fragments  are  elongated.  Studies  of  perforation  of  thin  metal 
plates  shows  dependence  of  the  ballistic  limits  on  the  impact  velocity  and 
the  mass/presented-arec  ratio.  Very  slender  projectiles  penetrate  and  per¬ 
forate  plates  more  readily  than  blunt  projectiles;  maximum  short  duration 

stresses  in  the  plate  are  dependent  on  the  impact  velocity  only.  If 

1/3  1/3 

damage  mechanisms  react  to  these  parameters  m  V  =  am  '  +  C  (or  V  =  con- 

stant)(See  Chapter  V,  page  86)  would  be  the  form  of  the  criterion.  The 
necessity  for  study  of  damage  by  slender  fragment  is  indicated  by  the 
conclusion  of  this  study  that  these  are  the  fragments  that  propagate  to 
the  larger  ranges. 
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